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PREFACE 


The  purpose  of  this  study  was  to  analyze  the  effects  of  mirror  mis¬ 
alignment  on  the  transverse  modes  and  beam  steering  of  an  unstable  laser 
resonator.  The  analysis  was  developed  for  any  general  unstable  resonator 
design  with  rectangular  apertures  and  did  not  allow  for  inclusion  of  a 
gain  medium.  The  final  result,  a  computer  program,  can  be  used  to  cal¬ 
culate  mode  eigenvalues  and  subsequently  the  intensity  and  phase  in  the 
plane  of  the  feedback  mirror  for  a  desired  mode.  The  slope  of  the  phase, 
from  which  a  beam  steering  angle  can  be  determined,  is  also  calculated 
for  the  lowest  loss  mode.  The  resulting  tilt  in  the  phase  front  is  due 
to  diffraction,  and  is  consequently  a  beam  steering  angle  additional  to 
the  geometric  misalignment.  The  code  is  basically  an  extension,  or  mod¬ 
ification,  to  a  previous  computer  model  developed  by  J.E.  Rowley,  and 
follows  similar  work  done  by  P.  Horwitz. 

Although  a  major  portion  of  this  work  may  be  found  elsewhere,  speci¬ 
fic  details  and  applications  throughout  the  text  are  generally  not  avail¬ 
able.  They  are  provided  here  in  order  to  present  a  complete  and  clear 
progression  of  the  analysis.  A  large  number  of  equations  and  derivations 
are  required  since  the  topic  is  analytical  in  nature  rather  than  experi¬ 
mental.  This  sometimes  leads  to  trivial  substitutions  and  algebraic 
steps  being  included  for  the  sake  of  continuity,  but  it  is  hoped  that 
these  are  minimal.  Physical  interpretations  and  definitions  are  included 
where  possible  for  better  understanding. 

I  would  like  to  express  my  gratitude  to  my  advisor,  Lt.  Col.  John 
Erkkila,  for  his  time,  patience,  guidance,  and  especially  his  enthusiasm 
which  inspired  me  continually  throughout. 
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ABSTRACT 


The  integral  equation  that  describes  mode  structure  of  an  unstable 
resonator  with  rectangular-  apertures  is  developed  from  scalar  diffraction 
theory.  This  equation,  modified  to  account  for  misalignments,  is  solved 
by  applying  the  asymptotic  methods  developed  by  Horwits.  A  second  order 
approximation  of  the  method  of  stationary  phase  is  then  employed  to  cal¬ 
culate  phase  and  intensity  values  for  all  points  in  the  output  plane. 

The  phase  front  is  also  curve  fitted  to  a  straight  line  over  the  geomet¬ 
rical  region  for  the  lowest  loss  mode.  From  the  slope  of  the  straight 
line,  a  direction  of  propagation  can  be  attributed  to  the  wave.  This  is 
a  diffracted  beam  steering  angle  and  is  additional  to  the  geometric  steer¬ 
ing  angle  (i.e.,  the  beam  steering  angle  due  to  the  geometric  misalign¬ 
ment  of  either  or  both  mirrors). 

Plots  of  intensity  and  phase  for  vaxious  degrees  of  misalignments 
are  presented  as  results  of  a  computer  program  that  utilizes  the  derived 
expressions.  Also  included  axe  graphs  of  the  phase  slope  versus  mirror 
misalignment. 


MODE  ANALioIS  IK  A 


MISALIGNED  UNSTABLE  RESONATOR 


I.  Introduction 


Background 

In  any  laser  cavity  there  are  two  inodes  of  oscillation:  longitudi¬ 
nal  (or  axial)  and  transverse  (or  radial).  The  longitudinal  modes  cou¬ 
pled  with  a  specific  transverse  mode  determine  the  frequency  of  radiation 
emitted,  while  the  transverse  modes  alone  determine  the  intensity  distri¬ 
bution  across  the  laser  beam.  In  either  case,  the  existing  modes  can  be 
found  by  applying  the  appropriate  boundary  conditions  to  the  wave  equa¬ 
tion.  In  a  stable  resonator  with  spherical  mirrors  and  edge  effects 
(diffraction  effects  due  to  the  finite  size  of  the  mirrors)  neglected, 
the  transverse  solutions  to  the  wave  equation  become  Hermite -Gaussian  or 
Laguerre-Gaussian  functions  for  mirrors  with  either  rectangular  or  cir¬ 
cular  symmetry  respectively  (Ref.  1:1324). 

For  the  unstable  resonator,  edge  effects  must  be  included  since  the 
output  is  a  diffraction-coupled  beam  passing  around,  rather  than  through 
the  output  mirror.  An  iterative  approach  to  finding  the  resonant  modes 
of  the  cavity  accounts  for  these  edge  effects  and  was  introduced  by  Fox 
and  Li  (Ref.  2).  This  method  is  to  first  assume  an  initial  field  distri¬ 
bution  in  the  resonator  and  then  apply  the  Fresnel-Kirchhoff  diffraction 
formula  to  the  problem  of  infinite  strip  mirrors. 

Figure  1-1  shows  differences  between  a  stable  resonator  and  an  un¬ 
stable  resonator.  Any  ray  originating  in  the  stable  resonator  and  strik¬ 
ing  one  of  the  mirrors  will  remain  inside  the  cavity,  even  for  an  infi- 


g.  =  1  -  £  i  =  1,2 
6i  «. 
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(l.D 


Here  L  =  cavity  length  and  =  radius  of  curvature  of  the  i  ‘  mirror. 
Also,  Ft,  is  defined  as  positive  for  concave  mirrors  and  negative  for 


0  < 


s  1 


(1.2) 


'Unstable  resonators  however,  are  characterized  by  either  g,g  >  1  or 

12 

g,g,  <  0  • 

The  equivalent  lens  train  of  a  typical  ’unstable  resonator  is  shown 
in  Figure  1-2  .  As  stated  earlier,  an  initial  field  distribution  is  as¬ 
sumed  at  P  and  propagated  through  one  round  trip  to  using  the 

Fresnel-Kirchhoff  diffraction  formula.  Kith  the  restriction  that  the 
field  reproduce  itself  after  the  round  trip  and  introducing  phase  lag  due 
to  mirror  curvature,  the  resulting  integral  equation  is 


l 


Here  x  is  the  spatial  coordinate,  y  is  a  dummy  variable  of  integra- 


tion,  i  =  &  ,  M  =  cavity  magnification,  g(x)  is  the  field  distribution 
on  the  output  mirror,  and  t  is  defined  in  equation  (B.10).  The  trans¬ 
verse  field  distribution  of  the  resonator  modes  are  then  given  by  the 
eigenfunctions  of  the  integral  equation.  The  multiplicative  constant,  v, 
is  the  eigenvalue  associated  with  the  eigenfunctions  and  gives  the  dif¬ 
fraction  loss  and  phase  shift  of  the  mode  (Ref.  1:1325). 

Rowley  (Ref.  3)*  following  the  analysis  of  Horwitz  (Ref.  4),  developed 
a  computer  code  to  numerically  solve  the  integral  equation.  He  assumed 
that  the  field  on  the  mirror  before  the  round  trip,  g(y),  consisted  of  a 
unit  amplitude  cylindrical  wave  plus  a  series  of  edge  diffracted  waves, 
and  is  given  by 

K 

g(y)  =  1  +  n|a  cnHn(y)  (1.4) 

(Ref.  3*17).  Therefore,  the  field  on  the  mirror  after  the  round  trip, 
g(x),  would  differ  only  by  the  multiplicative  constant  v  . 

Letting 


cnHn(y)  =  ^(y)  +  VnW 


(1.5) 


the  integral  equation  to  be  solved  is 


N 


v  Cl  +  n|1  C^Pn(x)  ♦  bnGn(x)]}  = 


S  / e"it(y '  -  JLiv»m  - 


dy  (1.6) 


-1 


The  method  of  stationary  phase  is  used  to  approximate  this  integral. 


It  states  that  an  integral  of  the  forn 


I 


(1.7) 


car.  be  expressed  as  a  series  when  t  is  large  and  q(y)  is  a  slowly  vary¬ 
ing  function.  The  first  two  terms  of  the  series  are  given  by 


T 


e-itp(y0) 


+ 


q(b)  -itp(b) 

pTb) 


q(a)  -itp(a) 
P'UJ 


(1.8) 


(Ref.  3;19;  14:1073).  This  is  a  first  order  approximation  to  the  inte¬ 
gral  equation  and  is  used  in  Chapter  4  for  determination  of  eigenvalues. 
A  second  order  approximation  is  also  given  in  Chapter  4  and  is  required 
for  intensity  and  phase  calculations  outside  the  mirror  edges.  In  equa¬ 
tion  (1.8)  the  prime  indicates  the  derivative  of  the  function  with  re¬ 
spect  to  y,  and  y0  is  the  stationary  phase  point,  i.e., 


p'(y0)  =  0 


(1.9) 


Before  equation  (1.6)  can  be  solved,  explicit  forms  for  Fn(y)  and 
cn(y)  axe  needed.  They  are 


Fn(y) 


exp[-it(l  -  y/Kn)/''Hn-1~1 
1  -  y/Mn 


(1.10) 


and 
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(1.11) 


G  (y)  =  -  /Hn-i  expf_-it(l  y/Kn)/^Mn_.,  ‘ 
n  V  4iTrt  1  +  y/M" 

(Ref.  3^20;  4; 1331) »  where  the  variable  M  is  defined  in  equation 
(3.12).  Horwitz  determined  the  functional  form  of  the  Fn's  and  CR's 
by  an  asymptotic  expansion  of  the  integral  equation  (Ref.  4).  The  spe¬ 
cific  set  of  functions  were  then  chosen  for  their  reproducibility,  i.e., 
after  making  the  stationary  phase  approximation  to  equation  (1.6)  the 
same  set  of  functions  are  arrived  at.  Equating  coefficients  then  leads 
to  a  2N  +  1  degree  polynomial  in  v,  which  is  solved  by  a  generalized 
root  finding  routine.  Once  the  eigenvalues  have  been  determined,  the 
field  is  calculated  by  first  specifying  a  particular  eigenvalue  or  mode 
and  then  calculating  the  constants  and  b^  (a  detailed  solution  is 

left  to  Chapter  4). 

A  limitation  to  this  analysis  is  the  fact  that  it  is  only  valid  for 
cavities  with  perfectly  aligned  mirrors.  Since  mirror  alignment  is  crit¬ 
ical  to  the  effectiveness  of  any  laser  system,  it  is  desirable  to  be  able 
to  determine  the  sensitivity  of  the  system  for  a  specified  misalignment. 
This  sensitivity  is  measured  in  terms  of  the  power  out,  beam  steering, 
and  mode  distortion  (i.e.,  the  beam  quality).  Several  investigations  in¬ 
to  mode  distortion  due  to  mirror  misalignment  have  been  reported  (Ref.  5> 
6,7,8),  however,  only  a  few  have  dealt  with  the  problem  of  beam  steering 
(Pef .  8). 

Objectives 

The  purpose  of  this  study  is  to  analyze  the  effects  of  mirror  mis¬ 
alignment  on  the  transverse  modes  and  beam  steering  of  an  unstable  laser 
resonator.  The  analysis  is  developed  for  a  bare  strip  resonator  with 
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rectangular  apertures  following  similar  work  done  by  Horwitz  (Ref.  5). 

A  computer  program  was  written  (a  modified  version  of  a  code  devel¬ 
oped  by  Rowley)  to  facilitate  the  study.  The  computer  model  calculates 
mode  eigenvalues  and  subsequently  evaluates  the  intensity  and  phase  in 
the  plane  of  the  feedback  mirror  for  a  desired  mode.  The  slope  of  the 
phase,  from  which  the  direction  of  propagation  can  be  determined,  is  al¬ 
so  calculated  for  the  lowest  loss  mode.  The  resulting  tilt  in  the  phase 
front  is  due  to  diffraction,  and  is  a  beam  steering  angle  additional  to 
the  geometric  misalignment. 


Assumptions 

The  basic  assumptions  of  Reference  3  2re  still  invoked  since  the 
analysis  is  an  extension  of  this  previous  work.  In  addition  to  these  is 
the  added  restriction  of  small  misalignments  (Assumption  5)-  They  are 
(Ref.  3=2)  : 

1.  That  scalar  diffraction  theory  can  be  used  if  a)  the  diffract¬ 
ing  aperture  is  large  compared  to  wavelength  and  b)  the  diffracted 
fields  are  not  observed  toe  close  to  the  aperture.  The  first  ccnstraint 
is  easily  attained  at  optical  wavelengths  and  the  second  is  achieved  at 
moderate  cavity  lengths. 

k.  That  the  diffraction  integrals  and  mode  eigenfunctions  are  sep¬ 
arable  allowing  a  1-D  strip  resonator  to  be  used  in  the  foregoing  anal¬ 
ysis.  When  the  mirror  separation  is  very  much  larger  than  the  mirror 
dimensions,  the  problem  of  the  rectangular  mirrors  reduces  to  a  one-di¬ 
mensional  problem  of  infinite  strip  mirrors  (Ref.  2:dj>7). 

3.  Diffraction  effects  from  only  the  feedback  mirror  are  considered 
if  the  larger  mirror  is  considered  infinite  in  comparison  with  the  beam 
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spot  size  on  that  mirror. 

9.  That  the  modes  in  the  strip  resonator  consist  of  a  fundamental 
cylindrical  wave  modified  by  a  finite  number  of  edge  diffraction  effects. 
This  assumption  is  supported  by  early  analysis  of  unstable  resonators 
(Ref.  9:279). 

5.  That  the  optical  axis  (the  line  joining  the  centers  of  curvature 
of  the  two  mirrors)  does  not  pass  too  close  (within  a  Fresnel  zone)  to 
the  edge  of  the  small  mirror  (Ref.  5:167). 

Procedure  and  Organization 

In  Chapter  2  the  Fresnel-Kirchhoff  diffraction  formula  is  modified 
to  include  phase  lag  due  to  mirror  curvature  for  a  strip  resonator;  and 
in  Chapter  3  the  resulting  integral  equation  is  generalized  to  include 
effects  of  small  misalignments.  Chapter  “4  then  deals  with  its  solution 
using  the  stationary  phase  approximation.  In  Chapter  5  the  diffracted 
beam  steering  angle  is  presented  and  the  magnitude  compared  to  the  geo¬ 
metric  beam  steering  angle,  while  the  results  and  conclusions  are  left 
to  Chapter  6. 


I I .  Development  of  the  Integral  Equation 
from  the 

Fresnel-Kirchhoff  Diffraction  Formula 


Equation  (1.3)  is  derived  by  first  reducing  the  problem  of  rectan¬ 
gular  mirrors  to  a  one-dimensional  problem  of  infinite  strip  mirrors. 
This  analysis  follows  the  work  of  Fox  and  Li  (Ref.  2) .  The  next  step  is 
to  account  for  phase  lag  due  to  mirror  curvature,  which  can  be  found  in 
References  3  and  4  .  Although  this  derivation  can  be  found  in  the  ref¬ 
erences  cited,  it  is  included  here  for  easy  access  and  completeness  of 
the  study. 


1-D  Analysis  (Ref.  2; 436-457,  464-436) 

In  a  Cartesian  coordinate  system,  the  field  Ep  due  to  an  illumi¬ 
nated  aperture  A  is  given  by  the  surface  integral 

Ep(x.y)  =  f Ea(x,y)  e~lkr  (1  +  cosO)  d5  (2.1) 

{ 

where  Ea  is  the  aperture  field,  k  is  the  propagation  constant,  r  is 
the  distance  from  a  point  on  the  aperture  to  the  point  of  observation, 
and  0  is  the  angle  which  r  makes  with  the  unit  normal  to  the  aper¬ 
ture.  From  Figure  2-1,  where  propagation  is  from  right  to  left,  this  can 
be  written  an 


E(xJ,yi)  = 


c  i 

*//' 


.ya) 


-ikr 


(1 


+  -) 
r 


dx„  dv_ 


r 


(2.2) 


Fig.  2-1.  Geometry  of  rectangular  plane  mirrors. 


Also  X  is  the  wavelength  of  the  radiation  and  r  is  the  distance  de¬ 
fined  by 


r  -  *  (x2  -  x,)'  -  (ys  -  (2-3) 

If  L  is  large  compared  to  the  dimensions  of  the  diffracting  aperture, 
mirror  2,  then  the  binomial  expansion  of  the  above  square  root  can  be 
approximated  by  retaining  only  the  first  two  terms  of  the  expansion,  i . e 

r  *  L[1  *  |(X4  -  *1)*  -  hi i  :  I  (2.4) 

where  L2  was  factored  out  of  the  square  root  before  the  expansion.  If 
again  L  is  considered  large  compared  to  dimensions  a  and  c,  then 
cos0  *  r-  ~  1  .  It  is  seen  here  that  even  though  r  ~  L,  the  exponential 


Fig.  2-2.  Typical  unstable  resonator  geometry. 


will  be  more  sensitive  to  this  round-off,  therefore  equation  (2.4)  is 
substituted  into  equation  (2.2)  to  become 

c  a 

=  1JT~'  /  f  E(xa*y8)  e 

~c  -a  (2.5) 

The  resonant  modes  of  any  laser  cavity  are  characterized  by  repro¬ 
ducible  fields  after  only  one  round  trip  through  the  resonator.  From 
Figure  2-2  this  is  given  by 


YE(x;,yp  =  E(xs, 


(2.6) 


Here  E(xg,y|J  is  the  field  after  the  round  trip  at  F  ,  E(xg,y,)  is  the 
original  field  at  P,  ,  and  Y  is  a  complex  constant.  When  the  mirrors 
have  rectangular  geometry  the  fields  are  best  represented  in  a  Cartesian 
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coordinate  system  where  the  two  orthogonal  components  become  separable. 


E(x,y)  =  U(x)  U(y) 


(2.7) 


Substitution  of  equation  (2.7)  into  equation  (2.5)  yields 


U(xi)U(yi) 


IK/ 

ZL/*2 


ik/  \  c 

2L'y2  ”  ‘Vi' 
e  dv„ 


"a  _c  (2.8) 
Restricting  the  analysis  to  a  one-dimensional  strip  resonator  with  the 
integral  defined  over  the  region  of  the  diffracting  aperture,  mirror  2, 
the  equation  for  the  one-way  diffraction  of  a  wave  from  F_  to  be¬ 

comes 


U(x 


ikL  f 

,> -7*.'—  / 

-a_ 


ik/  u 

U(x2)  e  dx2 


(2.9) 


Similarly,  the  one-way  diffraction  formula  for  propagation  from  Fj  to 


P'  is 


(2.10) 


Introduction  of  Phase  Lag 

For  equations  (2.9)  and  (2.10)  to  adequately  represent  propagation 
in  the  resonator,  the  phase  lag  due  to  mirror  curvature  must  be  accounted 
for.  Considering  equation  (2.9),  or  propagation  from  mirror  2  to  mirror 
1  in  Figure  2-2,  the  total  phase  lag  at  some  position  x  is 


0(x)  =  k[Al  -A(>1)  +  A_  -  a(x2)J 


(2.11) 


0(x)  -  k[R1  -  h*  -  x*  +  R2  -  -  x|  ] 


(2.12) 


The  maximum  free  space  region  between  plane  and  the  1  mirror  is 

A^,  while  A(x^)  is  the  free  space  region  between  plane  and  the  i u“ 

mirror  at  some  distance  Xy  .  It  should  be  emphasized  that  this  is  the 
phase  lag  for  propagation  in  one  direction  only. 

We  can  assume  the  paraxial  approximation  to  be  valid  over  the  entire 
mirror  if  the  radius  of  curvature  is  large  and  the  physical  dimensions 
are  small  compared  to  the  resonator  length.  Therefore,  if  the  above 
square  roots  are  expanded  in  a  binomial  series  as  before,  only  the  first 
two  terms  are  retained. 


h  -  Xj/R®  *  1  -  x^/2R‘ 


(2.13) 


Equation  (2.12)  is  now  written 


0(x)  =  k(x^/2R1  +  x|/2Re) 


(2.14) 


Writing  l/Rp  in  terms  of  the  g  parameter  from  equation  (l.l),  the 
phase  lag  takes  the  form 


0(x)  =  afxjd  -  Sl)  +  x|(i  -  g2)] 


(2.15) 


This  phase  term ,  seen  to  be  positive,  is  actually  a  phase  advance- 


13 


* 


ment;  but  when  added  to  the  negative  phase  of  equation  (2.9)  it  gains  the 
awkward  notation  of  phase  lag.  The  total  phase  for  the  one-way  diffrac¬ 
tion  of  a  wave  as  it  applies  to  a  resonator  with  spherical  mirrors  is 
then 


^total  =  '  +  x*  -  2x^  '  xi(1  -  61)  -  x"(1  '  6a)-' 

(2.16) 

After  combining  some  terms,  the  diffraction  formula  for  propagation  from 
to  Fx  is  given  by 


_  il-i  r 

\/5e’  2  J 

-a„ 


U<x>)  =  \ITl 

Similarly  equation  (2.10)  becomes 


ik/  £  ,  p  „ 

n,  v  -  2L(xig!  +  %  -  "iX2}  , 
U(x2)  e  dx 


s 

(2.1?) 


/a 

i -  -  f  -  ~(x'fc  +  y'  -  2x'x’) 

rr  2  I  ,N  2LIX1  X2  Ss  XiV  , 

U(X£}  =  JTL  *  J  U(xi}  6  dxi 

-a. 


(2.16) 

If  the  reproducibility  argument  of  equation  (2.6)  is  invoked,  i.e., 


Y' U(x ^)  =  U(xj) 


(2.19) 


then  equations  (2.17  -  2.19)  can  be  combined  to  yield  the  round  trip  dif¬ 
fraction  formula  as  it  applies  to  a  laser  resonator  with  spherical  mir¬ 
rors  and  rectangular  apertures.  In  the  equation  below,  the  constant 
phase  term  has  been  absorbed  into  the  complex  constant  y'  to  become  y  . 
The  integral  equation  is 
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a,  a? 

Y“Ui)  =  ±J  J  «u 


y  '  *f  y  1  'r  —  2x  1  Y  1  ) 
2L  i  6t  y‘£  ^xi  s' 


-a3  -a. 


-  M(X»S1  +  X^  -  2xlXa}  .  ,  , 

e  dxx  dxp 


(2.20) 


This  equation,  although  derived  from  the  geometry  of  Figure  2-2,  is 
equally  valid  for  the  resonator  of  Figure  1-2  and,  in  general,  any  laser 
resonator.  This  can  be  seen  from  the  sign  convention  chosen  for  the  ra¬ 
dius  cf  curvature,  Rj_  .  In  the  above  example  R  }  and  R_  are  defined 
as  positive  (refer  to  page  2)  leading  to  a  direct  substitution,  from 
equation  (2.14)  to  equation  (2.15),  of  l/R^  in  terms  of  gq  .  From 
Figure  1-2  however,  equation  (2.14)  takes  the  fora 


0(x)  =  k(x^/2Ra  -  x|/2Ea) 


(2.21) 


But  with  l/R  =  (g  -  l)/L  ,  the  phase  again  takes  the  form  of  equation 
(2.1 5).  Thus  equation  (2.20)  is  a  general  result  where  the  sign  cf 
has  been  absorbed  into  the  g  parameter. 

The  next  step  is  to  consider  diffraction  effects  from  the  feedback 
rirror  only.  This  permits  the  limits  of  integration  over  F  tc  go  to 
infinity  since  the  beam  spot  size  on  mirror  1  is  assumed  small  compared 
to  mirror  dimensions.  If  in  addition  to  setting  the  limits  of  integra¬ 
tion  over  mirror  1  to  i  05  ,  we  let 


x;  = 


(2.22) 
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the  integral  equation  becomes 


cL^  oo 


YU(x)  = 


iff 

-ap  -  •» 


U(y)  e 


-  ir(xi6i  +  x%  '  2x!x> 


-  ^(x‘6l  +  y  gj  -  2xjy) 

e  dx  a  dy 


(2.23) 


The  justification  for  equations  (2.22)  is  seen  by  following  a  single 
ray  from  mirror  2  to  mirror  1  and  then  back  again.  As  the  ray  strikes 
mirror  1  the  point  of  incidence  must  be  the  sane  as  the  point  of  reflec¬ 
tance,  i.e.,  xa  =  x'^  .  However,  the  ray  leaving  mirror  2  will  in  gener¬ 
al  be  displaced  from  the  original  position  giving  x£  /-  x^  .  Therefore, 
if  xj,  is  the  x  position  at  which  the  field  is  to  be  calculated  after 
the  round  trip  and  x2  =  y  (where  y  is  a  dummy  variable  of  integration 
over  the  diffracting  aperture),  equation  (2.20)  becomes  equation  (2.23) 
above. 

The  interior  integral  is  extracted  with  its  evaluation  relegated  to 
Appendix  A;  the  result  being  the  complete  kernel  of  the  integral  equation. 


ITT 


K  = 


ZXLgz 


ZXL; 


l(2§iS-  -  i)U2  +  yE)  -  2*y~1 


(2.2b-) 


Now  if 


(2.25a) 

(2.25  b) 

are  introduced,  the  round  trip  diffraction  formula  becomes 


G  =  2 -  1 

and 
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Y’J(x)  = 


(2.26) 


^_£(x2  +  y»)  -  2xyl 

az 


dy 


All  quantities  are  the  same  as  defined  previously  with  the  addition  of 
.  Here  Fg  is  u'ne  ordinary  Fresnel  number  of  the  smaller  mirror. 

The  ordinary  Fresnel  number  is  physically  interpreted  as  the  additional 
path  length  per  pass  in  half  wavelengths  for  a  ray  ^raveling  from  one 
mirror's  center  to  the  other  mirror's  edge,  compared  to  one  traveling 
from  mirror  center  to  mirror  center  (Ref.  Il:159~l6l). 

Normalizing  the  coordinate  system  such  that  a_  =  1  ,  the  resulting 
integral  equation  is 


yU(x)  =  v'lF  /  U(y) 


3-iTrFLg(xE  +  ys)  -  2xy]  dy 


(2.2?) 


Further  simplification  is  obtained  by  defining 


Heq  =  f(K  -  1/M)  (2.28) 

and 

U(x)  =  g(x)  (2.29) 


Several  different  interpretations  for  the  equivalent  Fresnel  number, 
Keq  ,  are  available  (Ref.  11:159-161;  12:360).  Here  the  equivalent 
Fresnel  number  is  defined  as  the  distance,  in  half  wavelengths,  between 
the  outer  edge  of  the  output  mirror  and  the  nearest  point  on  the  outgoing 
geometrical  wave  when  that  wave  just  touches  the  mirror  center  (Ref.  12: 
360).  That  geometrical  wave  is  assumed  to  be  a  cylindrical  wave  of  the 


1? 


form 


-inN_nx2 
e  ecl 


(2.30) 


Substitution  of  equations  (2.28)  and  (2.29)  into  equation  (2.2?)  yields 


-a 


-  iTtF[g(>:2  +  y2)  -  2xy] 

dv  (2.31) 


After  some  manipulation,  detailed  in  Appendix  B,  this  equation  sim¬ 


plifies  to  the  final  form. 


Ill .  The  Integral  Equation  Appropriate 
to  the 

Misaligned  Resonator 

In  Chapter  2  the  integral  equation  was  developed  for  the  case  of  a 
perfectly  aligned  resonator.  It  has  been  shown  (Ref.  y.  167-168;  10:2241- 
2242)  that  the  equation  appropriate  to  the  general  misaligned  resonator 
differs  from  the  usual  one  only  in  the  limits  of  integration.  This  is 
true  regardless  of  which  mirror  is  tilted. 

The  limits  of  integration  are  determined  by  the  angle  through  which 
the  mirror  is  tilted,  which  mirror  is  misaligned,  and  the  cavity  geome¬ 
try.  The  results  of  this  section  (i.e.,  the  limits  of  integration)  per¬ 
tain  only  to  the  geometry  of  Figure  and  cannot  be  generalized  to  in¬ 
clude  other  configurations. 

Effects  of  Mirror  Misalignment 

This  analysis  is  basically  a  geometry  problem  with  the  development 
restricted  to  misalignment  of  the  feedback  mirror. 

In  Figure  3~1  the  optic  axis  of  the  perfectly  aligned  resonator  is 
the  line  bcde.  When  mirror  2  is  tilted  around  point  c  by  an  angle  0,  the 
optic  axis  of  the  resonator  becomes  hgfe.  The  center  of  curvature  of 
mirror  1  is  at  e,  while  the  centers  of  curvature  of  mirror  2  before  and 
after  the  tilt  are  d  and  f,  respectively.  The  angle  between  the  new  op¬ 
tic  axis  and  the  old  is  o>.  For  small  angles  0,  oj  is  given  by 


Mirror  1 


Fig.  3-1.  Geometry  of  am  unstable  laser  resonator. 

This,  however,  neglects  the  fact  that  R£  by  convention  is  negative. 
Therefore,  substitution  of  a  negative  E£  and  rewriting  in  terms  of  the 
g  parameter,  equation  (3.1)  becomes 


CD  =  G 


61  ~  1 

1  ~ 


(3.2) 


The  degree  of  tilt,  as  expressed  in  distance  across  the  feedback 
mirror,  is  given  by  the  line  eg”  .  This  distance  is 


eg  -  -  L* )  «  a)(R,  -  L)  (3-3) 

If  we  again  write  R3  in  terms  of  g3  and  use  the  results  of  equation 
(3.2) ,  it  becomes 


2  ( 


The  problem  now  is  to  express  the  distance  eg  in  terms  of  the  mirror 


dimension  a 

2 


eg  _  LG  Sj 

a2  a2  6i6a  -1 


Prom  equation  (2.2_5)  this  is  ’written 


eg  =  J.  ap/x 
a2  2F  glg;  -  1 


With  the  following  definitions 


(3.5) 


(3.6) 


N'eq  =  |  F(M  "  V?0 

and 

(g262)z  +  (StSg  ~  O2 
(gxg2)z  ~  (gjgg  -  !)z 


(3.7a) 


(3.7b) 


equation  (3*6)  takes  the  form 


=  Cg  _ 


a, 

r 


e  m  +  l 


N 


eq 


M  -  1 


(3.8) 


When  the  angle  6  is  small,  the  tilted  resonator  is  equivalent  to  an 
aligned  resonator  with  an  asymmetric  mirror  relative  to  the  resonator 
axis  (Ref.  10:2242),  refer  to  Figure  .  Since  the  origin  of  the  x- 
axis  is  defined  as  the  point  of  intersection  of  the  optic  axis  and  mir- 


ror  2,  the  extent  of  the  asymmetric  mirror  is  from  -1+6  to  1  +  6  . 
The  limits  of  integration  are  now 


a  =  -1  +  6 


(3-9a) 


and 


3  =  i+6 


(3-9b) 


Although  6  is  in  general  a  linear  function  of  0  for  any  cavity  con¬ 
figuration,  its  exact  form  must  be  determined  from  the  specific  geometry 
The  generalized  integral  equation  is  now  written  as 


vg(x) 


g(y) 


-it(y  -  x/M) : 


dy 


(3.lo) 


The  series  functions  Fn  and  Gn  also  have  minor  changes  due  to 
the  misalignment,  i.e., 


and 


Fn(x)  =  - 


3  -  x/M 


M„_.  exp[_-it(3  -  x/Mn)7Hn_1] 

n 


(3.Ua) 


Gn(x)  = 


exp_-it(a  -  x/Mn) '/Mn_ 3 " 
a  -  x/Mn 


(3.Hb) 


(Ref.  5; 168-169) ,  where  again  M  is  the  cavity  magnification  and 


M 


n 


n 

L. 

k=0 


M-2k 

M 


(3.12) 


As  previously  stated,  6  is  a  function  of  which  mirror  is  tilted. 


If  mirror  1  had  been  misaligned  by  the  same  angle  9,  then  it  would  have 
been  found  (Ref.  8:579)  that 


eg 


6R  R 
i  s 


9L 


L  -  Rj  -  Re 


Sjgg  -  i 


(3.13) 


It  is  evident  from  equations  (3.4)  and  (3*13) »  that  the  geometry  of  Fig¬ 
ure  3-1  is  less  sensitive  to  misalignment  (by  a  factor  of  g1)  for  a  tilt 
of  mirror  2  vs.  the  same  tilt  of  mirror  1  . 


IV.  Solution  of  the  Generalized 


Integral  Equation 


In  Chapter  3  the  integral  equation  was  modified  to  include  effects 
of  small  misalignments.  This  chapter  is  concerned  witn  its  solution  by 
applying  the  method  of  stationary  phase. 

Since  the  problem  is  basically  algebraic  and  requires  numerous  man¬ 
ipulations,  the  scope  of  this  section  is  limited  to  a  detailed  outline 
of  the  solution. 

Stationary  Phase  Approximation 

The  integral  equation  to  be  solved  is 

V  6(x)  =  y?  J  s(,)  e'1^  -  dy  (*.1) 

a 

where  g(x)  is  the  field  distribution  on  the  output  mirror  after  the  round 
trip,  g(y)  is  the  original  field  distribution,  and  v  is  the  eigenvalue 
associated  with  the  eigenfunctions. 

Following  Chapter  1, 

N 

g(y)  -  1  +  l  cnH  (y)  (4.2) 

n-i 

where  it  was  assumed  that  g(y)  consisted  of  a  unit  amplitude  cylindrical 
wave  plus  a  series  of  edge  diffracted  waves.  The  basis  for  this  assump¬ 
tion  is  that  the  original  field  is  composed  of  the  primary  cylindrical 
wave  plus  diffraction  effects  from  the  previous  N  reflections.  Tne 
scries  terminates  with  Hj.(y),  the  last  function  to  have  an  effect  on  the 
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field.  The  addition  of  one  more  function  (or  reflection)  to  the  series 
would  add  only  to  the  amplitude  and  not  the  spatial  distribution  or  shape 
of  the  field,  i.e.,  H^+^(y)  is  a  constant.  A  good  approximation  is  to 
let 


j.}  >  In  (25O  Kgq) 
In  M 


(4.3) 


(Ref.  4:1533). 

Combining  equations  (1.5),  (4.1),  and  (4.2)  yields 


K 

v  [l  +  ^  LanFn(x)  +  cnGn(x)]]  = 


Lar,Fn(y)  +  trGn(y)nidy 


(4.4) 


This  is  equivalent  to  equation  (1.6)  with  the  limits  of  integration 
changed  to  account  for  the  misalignment.  Substitution  of  the  Fn's  and 
Gn's  from  equation  (3.H)  allows  the  stationary  phase  approximation  to  be 
applied,  refer  to  equation  (1.8).  Defining  the  quantities 


*o  = 


e-it(y  -  x/M )2  dv 


(4.5a) 


a 


:n  =  \[¥  /  LanFn(y)  +  bnGn(y)]  e'iL(y  “  x/M)  dy  (4.5b) 


equation  <’4. 4)  becomes 


N  U 

v  (1  +  r|,  anFn<x)  +  bnGn(x)  Tn 


(4.6) 
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The  problem  in  to  now  apply  the  stationary  phase  approximation  to  the 
individual  ln's  and  then  sum  the  results.  Starting  with  IQ  and  com¬ 
paring  with  equation  (1.7),  it  is  seen  that 


and 


q(y)  = 

(4.7a) 

p(y)  =  (y  -  x/m)' 

(4.7b) 

Also,  from  the  stationary  phase  point  definition,  equation  ( 1 . 9) • 


yQ  =  x/m 


(*.8) 


Employing  equation  (1.6)  and  with  some  manipulation,  it  can  be  shown  that 


=  1  - 

V  4iirt 


-it(P  -  x/K)c  -it((X  -  x/M )  ’ 


fs  -  x/M 


a  -  x/M 


(4.9) 


When  written  in  terms  of  Fn  and  Gn  from  equation  (3*11),  this  becomes 


10  =  1  +  Fx(x)  +  Gj(x) 


(4.10) 


Applying  the  same  technique  to  Iir  where 


= 


v  I  +  b^^y);]  e 


-it(y  -  x/M)' 


dy 


(4.U) 


the  approximation  yields 
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Ij  =  a1F„(x)  +  \Gs(x)  +  1 1(x)La1FJ(t)  +  b^O)] 


+  GjfxJLajF^a)  +  bjG^a)] 


(4.12) 


One  more  iteration  is  necessary  in  order  to  show  a  general  trend.  Again, 
using  the  stationary  phase  approximation  to  evaluate  1^,  the  result  is 


I.  =  a/3(x)  +  bsG3(x)  +  F1(x)[a,F,(  3)  +  b^fc')] 


4  G1(x)La,F2(a)  +  b2G£(a)] 


(4.13) 


From  equations  (4.10),  (4.12),  and  (4.13)  the  resulting  summation  can  be 
generalized  to  include  all  In's,  i.e., 

N  K 

n^0  Jn  =  1  +  Fz(>-)  +  Vx)  +  CanFn+1(x)  +  V^x)] 


+  Fl(x)  ^  [^(p)  +bhGn(&)D 


Gj(x)  Z  CanFn(a)  +  bnGn(a)] 


(4.14) 


The  equivalence  of  equations  (4.6)  and  (4.14)  gives  rise  to 


v  [1  +  ^  [a^F^x)  +  bnGn(x)]}  =  1  +  F1(x)  +  Ga(x) 


riir,  i-Vn+jW  4  b,:Gn+i<x)] 


4  £V„<e>  4  V=n<B>l  4  G,<*>  J,  4  b„G„(a)] 

(4.1ji) 
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The  Polynomial  Equation 

From  equation  (4.15)  a  polynomial  equation  with  determinable  coef¬ 
ficients  is  developed.  This  eigenvalue  polynomial  (where  v  is  the 
eigenvalue)  is  of  order  2N  +  1  and  is  solved  numerically.  Each  root 
is  then  associated  with  a  different  resonant  mode  of  the  cavity. 

Referring  to  equation  (4.15)*  the  first  step  is  to  equate  coeffi¬ 
cients.  For  n  f  1  we  have 


Si  +  1  =  an  ■  V 


(4. 16a) 


vb  .  =  b  =  b.Tv 

n+i  n  n 


(4. 16b) 


Here  the  right  equality  is  obtained  from 


fn  -  fi 

Si+I  V  n 


(4.17) 


and  letting  n  =  N  -  1  ,  i.e., 


aa  V 


^-1  N-l  =  ‘V 


(4.18) 


The  next  step  is  to  equate  constant  terms,  which  gives 


v  =  1  +  Sfl^+l  +  V*N+1 


(4.  HO 


Here  the  x  dependence  has  been  dropped  since,  by  definition,  the  last 
functions  to  contribute  to  the  spatial  distribution  of  the  field  are 
Fn(x)  and  GN(x),  or  that  F^+j(x)  and  Gj.+^(x)  are  constant  for  all  x  . 
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Equating  coefficients  of  F1(x)  and  G.(x)  allows  us  to  write 


1  + 


N 

nL  CanF„(«  * 


(,4.20a) 


and 


vb 


1 


1  + 


K 

I  Ca,rn(a)  +  b  G  (a) j 
n=:  ‘  *•  n 


(4.20b) 


Substitution  for  an,  bn>  a  ,  and  fc1  from  equation  (4.16)  and  rewriting 
equation  (4.19)  we  obtain 


V'”  =  1  +  IV^-H  C^O)  +  ^Gn(P)]  (4.21a) 

bt.vN  =  l  +  Iv;;'n  [ajjFn(a)  +  bj;Gn(a)]  (4.21b) 

V  1  +  alI^N+l  +  \tCN+1  ,4.21c) 


In  the  above  equations  the  summations  are  understood  to  be  from  n  =  1 
to  n  =  N  .  Although  not  readily  apparent,  equations  (4.21)  can  be  com¬ 
bined  to  yield  a  polynomial  equation  in  v  of  order  2N  +  1  ,  The  ex¬ 
pression  is  left  to  Appendix  C,  where  it  is  seen  that  the  coefficients 
are  calculated  from  a  knowledge  of  Fn(a),  -rn  O ) ,  Gn(cr),  G  (  p)  ,  Fj,+1,  and 
Gj.+^.  Therefore,  by  specifying  the  quantities  M,  6,  and  N  (refer  to 
equations  (3*9)»  (3*H)i  (3>12),  and  (4.3)  ),  the  roots  of  the  polynomial 
are  eventually  determined. 

Once  the  eigenvalues  have  been  calculated,  the  constants  a  and 

n 

b^  are  determined  for  a  particular  mode  (again  refer  to  Appendix  C). 

The  resulting  field  is  then  evaluated  at  incremental  positions  across  the 
feedback  mirror,  i.e., 
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4.21) 


N 

g(x)  =  1  +  L  La  F  (x)  +  b  Z  (>:) 

n=i  n  r‘  n  r- 

Given  the  field,  g(x),  the  phase  is  easily  calculated  fron 

jZf(x)  =  arctan(^)  (4.23) 

where  g(x)  is  the  complex  number  x  +  iy  . 

Second  Order  Approximation 

A  higher  order  approximation  to  the  method  of  stationary  phase  is 
required,  due  to  the  singularities  involved,  whenever  y  approaches  the 
endpoints  of  the  integral.  For  IQ  (refer  to  equation  (4.8)  )  and  the 
case  of  a  perfectly  aligned  resonator,  this  occurs  when  x  approaches 
the  shadow  boundaries. 

Since  the  v's  were  determined  from  a  valid  first  order  approxima¬ 
tion  over  the  region  G  <  x  <  f3  ,  then  the  eigenvalues  are  acceptable  for 
all  x  .  This  is  justified  by  noting  that  v  is  a  constant.  The  series 

constants  a  and  b  are  also  valid  from  the  first  approximation,  re- 
n  n 

fer  to  Appendix  C,  where  it  is  seen  that  they  are  a  function  of  the  con- 
N — n 

stants  v  ,  F  (a),  Gn(a)»  FK+1,  and  •  This  implies  that  in  order 

to  calculate  the  phase  and  intensity  for  all  x  values,  a  new  approxima¬ 
tion  to  the  integral 

(5 

I  =  J q(y)  e-ltpl'-V)  dy  (4.24) 

a 

is  all  that  is  required. 
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Tne  higher  order  approximation  j  equation  (4.24)  can  oe  simplified 
by  letting 


Depending  on  the  location  of  the  stationary  phase  point,  y  ,  the  approxi¬ 
mation  takes  one  of  three  forms.  They  are:  1)  for  yQ  £  a 


I  =  u(p)  E*i 


•(f)]  - 


’Lv(a)] 


26) 


and  3)  for  yQ  ^  6 

I  =  u(£)  j^E^CvO)]  +  — j  -  u(a)  |^E*[v(a)]  +  ■ ~  j  (4.28) 

(Ref.  3: 29-30) .  In  the  above  equations,  E  is  the  complex  conjugate  of 
the  Fresnel  integral. 

Since  the  mathematics  are  again  quite  involved  and  tedious,  the  ap¬ 
plication  of  equations  (4.26-4.28)  will  not  be  given.  It  is,  however. 
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sufficient  to  say  that  the  method  of  solution  is  similar  to  that  of  the 
first  order  approximation. 
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V.  BEAM  STEERING 


The  major  effects  of  mirror  misalignment  in  unstable  resonators  are 
beam  steering  and  mode  distortion.  Mode  distortion  is  easily  verified 
by  observation  of  the  results  presented  in  the  next  chapter.  The  geomet¬ 
rical  beam  steering  angle,  Y  >  is  however  not  determined  due  to  the  math- 

g 

ematical  construct  of  the  misaligned  resonator,  i.e.',  the  misaligned  re¬ 
sonator  is  modeled  as  an  aligned  asymmetric  resonator  (see  Figure  5-l) . 
Additional  information  is  required  in  order  to  calculate  y  (refer  to 
Chapter  3)  since  this  analysis  depends  only  on  a  knowledge  of  6,  M,  and 

N  . 
eq 

A  diffracted  beam  steering  angle,  Y is  observed  by  comparison  of 
Figures  6-2,  6-4,  and  6-6.  The  phase  fronts  are  seen  to  shift  slightly 
with  variations  in  the  parameter  6.  Since  the  normal  to  the  phase  front 
determines  the  direction  of  propagation,  a  straight  line  curve  fit  of  the 
phase  is  desired. 

This  chapter  calculates  order  of  magnitude  quantities  for  Y  and. 

o 

9  (the  mirror  tilt  angle).  Also,  a  least  squares  curve  fit  is  discussed 
for  determination  of  Y^* 


Geometrical  Beam  Steering 

The  beam  steering  angle  due  to  the  geometric  misalignment  is  the 
angle  which  the  new  optic  axis  makes  with  the  old.  From  Figure  3-l 


Y 


g 


CD 


9 


T 


-  1 

- 


(5.D 


Combining  equations  (3-8),  (3-7),  and  (2.23)  ihe  tilt  angle  of  the  mirror 
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can  be  written  an 


e 


a. 

Rm  -  l)2 

4Lgx 

K  J 

6 


(5-2) 


Substitution  of  equation  (5*2)  into  (5»l)  gives 


in  terms  of  6. 


Y 


S 


s-.  -  1“ 

—  — 

nfp 

*(H  -  1)£' 

L1"  SiS2_ 

4bgl 

R 

6 


(50) 


In  order  to  make  a  comparative  analysis  later,  we  specify  M  =  2.0, 
9.6,  and  choose  L  to  be  2.0  meters.  From  equations  (B.6.)  and 
(2.2 5) 


gaga  =  i(M  +  1/M  +  2)  =  1.125  (5.4) 

Calculating  R^  in  terms  of  Ra  is  possible  by  substitution  of  equation” 
(l.l)  in  the  above  expression.  The  cavity  configuration  of  Figure  5“1  is 
also  desired  (i.e.,  Ra  is  positive  and  Rg  is  negative),  which  gives 

i  qrshryr  (5-5) 

Having  already  specified  R  an  positive,  requires 

(eLi^Rg  <  10  (5.6) 

or  Ra  <  16  meters.  Therefore,  letting  Rs  =  10  meters  the  value  of  Ra 

is  set  at  32  meters.  The  maximum  geometrical  beam  steering  angle,  y  , 

&max 

occurs  for  6  =  1.0  .  For  6  >  1  there  is  no  longer  an  axis  within  the 
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Fig.  5-1.  The  equivalent  asymmetric  cavity  of  the  misaligned  re¬ 
sonator  as  discussed  in  Chapter  3* 

resonator  with  the  symmetry  of  the  original  resonator,  and  there  is  no 
possibility  of  exciting  modes  characteristic  of  the  aligned  unstable  re¬ 
sonator  (Ref.  8:579-580).  With  a„  =  0.015  meters,  Ygmax  is  equal  to 
0.5  mrajl.  This  is  equivalent  to  a  mirror  tilt  of  approximately  1  mrad 
as  calculated  from  equation  (5*2). 

Beam  Steering  Angles  Due  to  Diffraction 

The  diffracted  beam  steering  angle  is  first  determined  by  finding 
the  slope  of  the  phase  over  the  geometrical  region.  Refering  to  Figure 
5-1,  the  region  is  seen  to  be  from  -M+6M  to  M+&M.  This  is  easily  ob¬ 
tained  by  noting  that  the  cavity  magnification  is  a  constant  depending 
only  on  Rlt  Rp,  and  L.  The  transverse  magnification  of  the  perfectly 
aligned  resonator  is  defined  as 
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where  the  dimensions  are  normalized  such  that  a  =  1.0  .  When  mirror  2 

2 

is  misaligned  the  equivalent  resonator  of  Figure  5~1  has  the  same  cavity 
magnification,  and  is  given  try 

n  -  rtV  (5.8) 

2 

or  x2  =  M  +  6M.  The  extent  of  the  geometrical  region  in  the  negative 
direction  can  similarly  be  shown  to  be  -M  +  6m. 

A  straight  line  curve  fit  of  the  phase  is  achieved  by  the  method  of 
least  squares.  Given  n  sets  of  points  (x,y),  where  y  is  the  phase 
in  radians  and  x  is  the  normalized  distance,  the  best  straight  line 
y  =  mx  +  b  is  determined  by  solving  the  two  normal  equations 


nn  +  m  I  xj  =  2  yj 
b  I  Xj  +  m  2  x^  =  L  Xj  y. 


(5-9) 


(Ref.  15:683),  where  the  summation  is  from  j  =  1  to  j  =  n  .  The  dif¬ 
fracted  beam  steering  angle  is  then 


tan  vd  «  Vd  =  ~ 


(5-10) 


The  quantities  d  and  r  are  the  change  in  y  and  x  respectively  in 
MKS  units,  i.e.,  d  is  obtained  from  e^^  =  e^  and  r  =  a2x.  Therefore, 
the  diffracted  beam  steering  angle  is 


Prom  the  results  of  the  next  chapter,  the  slope  (y/x)  is  determined 

from  the  phase  calculations  for  various  values  of  6  and  plotted  in 

Figures  6-7  and  6-8.  Analyzing  the  case  of  K  =  2  and  N  =9.6  (Figure 

eq. 

6-7) ,  (y/x)max  is  seen  to  be  approximately  0.14  .  Also,  from  equations 
(2.28)  and  (2.25) 


X  =  a„(M  -  1/M) 
a2  4gllNeqL 


(5.12) 


which  gives 


a„(M  -  1/M) 

brd^qLgj 


(5-13) 


Using  the  values  from  the  previous  cavity,  a£  =  0.015  m  and  L  =  2.0  m, 

y.  «  7  p.rad.  The  maximum  value  of  0.1 4  is  seen  to  occur  at  6  =  .225, 
amax 

therefore,-  y  is  equal  to  0.225  y  or  «  0.1  mrad.  It  is  clearly 
6  §max 

seen  that  the  diffracted  beam  steering  angle  is  approximately  7%  of 


the  geometrical  beam  steering  angle,  and  is  a  significant  contribution 
when  considering  propagation  over  long  distances. 


VI .  Results,  Conclusions 


and 

Recommendations 


Results 

The  main  result  of  this  study  is  the  development  of  the  computer 
program,  BARC  2,  as  listed  in  Appendix  E  along  with  a  description  of  in¬ 
put  variables.  This  code  calculates  intensity  and  phase  in  the  plane  of 

the  feedback  mirror  by  specifying  &,  M,  and  N  .  Also,  the  slope  (y/x  ) 

eq 

is  determined  for  the  lowest  loss  mode  only,  where  the  diffracted  beam 

steering  angle  is  related  to  the  slope  by  equation  (5.13). 

Plots  of  intensity  and  phase  for  the  first  three  modes  of  a  cavity 

with  M  =  2.0  and  N  =9.6  are  at  the  end  of  this  section  and  in  Ap- 

eq  r 

pendix  D.  For  various  6,  the  intensity  and  phase  distortion  is  clearly 
evident.  The  phase  of  the  lowest  loss  mode  (Figures  6-2,  6-4,  and  6-6) 
is  seen  tc  remain  relatively  uniform,  but  has  a  definite  slope  or  direc¬ 
tionality  associated  with  it  an  &  changes.  This  slope  is  plotted  in 
Figure  6-7  for  a  range  of  6  between  0.0  and  0.25  • 

For  the  cavity  configuration  chosen  in  Chapter  5.  the  diffracted 
beam  steering  angle  is  a  significant  portion  of  the  total  beam  steering 
angle  for  specific  values  of  6.  Again  referring  to  Figure  6-7,  for  6  = 

0.0H25  the  slope  is  0.1304  .  This  gives  y  «  6.5  prad  and  y  = 

o  S 

0.01125  y  5.6  p rad,  which  shows  that  y,  is  of  the  order  of  Y 

Cmax  o  g 

when  6  ;•  0.01  .  This  value  is,  however,  totally  dependent  on  cavity 
geometry  and  must  be  evaluated  for  each  resonator  as  discussed  in  Chap¬ 
ters  3  and  5-  More  importantly,  \  is  seen  to  vary  with  no  apparent 
regularity  as  6  varies. 
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Conclusions 


By  comparison  of  the  phase  and  intensity  plots  with  those  in  Refer¬ 
ence  5*  *he  basic  conclusion  reached  is  that  program  BARC2  produces  valid 
results.  Also,  for  the  case  of  a  perfectly  aligned  resonator,  i.e.  6  = 

0,  the  results  are  consistent  with  those  of  the  previous  study  (Ref.  3)« 

As  shown  in  Chapter  3»  the  geometrical  beam  steering  angle  is  a  lin¬ 
ear  function  of  the  mirror  tilt  angle  for  small  misalignments.  This  has 
been  confirmed  from  a  previous  analysis  by  Krupke  and  Sooy  (Ref.  8).  Al¬ 
so,  the  diffracted  beam  steering  angle  is  of  the  order  of  the  geometrical 
beam  steering  angle  for  6  «  1.0;  however  the  exact  limit  depends  on 
which  mirror  is  misaligned  and  the  cavity  geometry. 

The  irregularities  associated  with  the  diffracted  beam  steering 
angle  versus  mirror  misalignment  (Figures  6-7  and  6-8)  is  a  result  of  the 
structure  in  the  phase.  Had  the  feedback  mirror  been  illuminated  by  a 
plane  wave  of  infinite  extent,  the  diffracted  beam  steering  angle  would 
have  been  zero  for  all  6  since  the  mathematical  model  is  simply  a  trans¬ 
lation  of  the  diffracting  obstacle  (see  Figure  5~l).  Referring  to  Figure 
6-2,  when  the  mirror  is  translated  to  the  right  (the  dashed  lines  indi¬ 
cate  the- position  of  the  mirror)  it  intercepts  a  further  advanced  wave 
on  the  right  than  on  the  left,  and  a  definite  slope  is  observed  in  the 
phase  front.  Specifically,  at  6  ~  0.2  (again  refer  to  Figure  6-2)  the 
extent  of  the  mirror  becomes  -0.8  to  1.2  where  a  peak  in  the  phase 
is  intercepted  on  the  right  and  a  minimum  on  the  left.  This  produces  the 
maximum  beam  steering  angle  of  Figure  6-7. 

Further  inspection  of  Figure  6-2  would  require  Yj  bo  occur  at 

amax 

6  w  0.9  .  This  is  exactly  the  case  and  the  slope  at  this  point  is  0.519. 
The  reason  fcr  terminating  the  plots  in  Figures  6-7  and  6-8  at  6  =  O.25 
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is  due  to  the  excessive  amount  of  computer  time  required  to  generate  them. 
Therefore,  relative  values  for  the  diffraction  beam  steering  -angle  can  be 
estimated  from  the  phase  of  the  perfectly  aligned  rescnator. 


Recommendations 

The  computer  model,  developed  from  the  previous  analysis,  calculates 
intensity  and  phase  in  the  plane  of  the  feedback  mirror  for  the  case  of  a 
bare  strip  resonator.  A  beam  steering  angle  due  to  diffraction  is  also 
determined  for  the  lowest  loss  mode. 

An  obvious  extension  would  be  to  determine  beam  steering  angles  for 
all  possible  modes  of  the  cavity.  In  addition,  the  model  might  be  modi¬ 
fied  to  account  for  the  presence  of  a  saturated  or  non-uniform  gain  medi¬ 
um.  Since  the  analysis  was  restricted  to  mirrors  with  rectangular  aper¬ 
tures,  the  case  of  circular  mirrors  would  be  another  topic  for  investiga¬ 
tion  . 

Still  another  area  of  consideration  would  be  the  analysis  of  mirror 
tilt  and  its  effect  on  higher  order  aberrations;  or  the  change  in  the 
phase  slope  vs.  mirror  misalignment  curves  for  higher  Fresnel  numbers. 
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Fig.  6-3.  Intensity  plot  for  the  lowest  loss  mode  with  .6  = 


0.2  . 


Fig.  6-4.  Phase  plot  for  the  lowest  loss  mode  with 


42 


BIBLIOGRAPHY 


1.  Kogelnik,  H.  and  T.  Li.  "Laser  Beams  and  Resonators,"  Proceedings 

IEEE,  Vol.  54,  No.  10:  1312-1329  (Oct.  1966) . 

2.  Fox,  A.G.  and  T.  Li.  "Resonant  Modes  in  a  Maser  Interferometer," 

Bell  System  Technical  Journal,  Vol.  40,  No.  2:  453-488  (March  1961). 

3.  Rowley,  James  E.  Computer  Analysis  of  Modes  in  an  Unstable  Strip 
Laser  Resonator.  MS  Thesis.  V’right-Pattersor.  AFB,  Ohio:  School  of 
Engineering,  Air  Force  Institute  of  Technology,  Dec.  I98O. 

4.  Horwitz,  P.  "Asymptotic  Theory  of  Unstable  Resonator  Modes,"  Jour¬ 
nal  of  the  Optical  Society  of  America.  Vol.  63.  No.  12:  I528-I5H-3 

(Dec.  1973). 

5.  - .  "Modes  in  Misaligned  Unstable  Resonators,"  Applied  Optics. 

Vol.  15,  No. 1:  167-178  (Jan.  1976). 

6.  Perkins,  J.F.  and  C.  Cason.  "Effects  of  Small  Misalignments  in  Empty 

Unstable  Resonators,"  Applied  Physics  Letters.  Vol.  31  *  No.  '}■  198- 

200  (Aug.  1977). 

7.  Hauck,  R. ,  H.P.  Kortz,  and  H.  Weber.  "Misalignment  Sensitivity  of 

Optical  Resonators,"  Applied  Optics.  Vol.  19,  No.  4:  598-601  (Feb. 

1990). 

8.  Krupke,  W.F.  and  W.R.  Sooy.  "Properties  of  an  Unstable  Confocal  Re¬ 

sonator  C0„  Laser  System,"  IEEE  Journal  of  Quantum  Electronics. 
Vol.  QE-5,  No.  12:  575-586  (Dec.  1969). 

9.  Siegman,  A.E.  "Unstable  Optical  Resonators  for  Laser  Applications," 

Proceedings  of  the  IEEE.  Vol.  53»  No.  3:  277-287  (March  1965)- 

10.  Sanderson,  R.L.  and  W.  Streifer.  "Laser  Resonators  with  Tilted  Re¬ 
flectors,"  Applied  Optics,  Vol.  8,  No.  11:  2241-2248  (Nov.  1969) • 

11.  Siegman,  A.E.  and  R.  Arrathoon.  "Modes  in  Unstable  Optical  Resona¬ 

tors  and  Lens  Waveguides,"  IEEE  Journal  of  Quantum  Electronics,  Vol. 
QE-3,  No.  4:  I56-I63  (April  I967)* 

12.  Siegman,  A.E.  "Unstable  Optical  Resonators,"  Applied  Optics,  Vol. 

13,  No.  2:  353-367  (Feb.  1974). 

13.  Beyer,  W.H.  CRC  Standard  Mathematical  Tables  (25th  Edition).  Boca 
Raton,  Florida:  CRC  Press,  Inc.,  1980. 

14.  Butts,  R.R.  and  P.V.  Avizonis.  "Asymptotic  Analysis  of  Unstable  La¬ 

ser  Resonators  with  Circular  Mirrors,"  Journal  of  the  Optical  Soci¬ 
ety  of  America.  Vol.  68,  No.  8:  1072-1076  (Aug.  1978) . 


46 


APPENDIX  A 


Appendix  A  is  devoted  to  the  evaluation  of  the  definite  integral  in 
equation  (2.23). 


OC 

ij 


ik/  2 


Eg£  -  2xyx)  -  |£(xfga  +  y2g, 
e 


-  2xay) 


dx  ( A .  1 ) 


The  exponentials  can  be  combined  to  yield 


ik  ( 
~  2L  Bs^ 

XL  e 


»s  +  »a)  f -i 


-  |fC2x2g;i  -  2xa(x  +  y)] 


dx. 


(A. 2) 


Completing  the  square  in  the  exponent,  the  integral  becomes 


vu 

.  .  (xa  +  iktxnil3  f  .  it[x=E  -  y  (x  +  v)  +  4  4 
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(A. 3) 


Instead  of  ( A . 3) »  we  write  the  integral  as 


i  r  ikr 

XL  exp[”  2L-6^ 


,(x2  +  y2) 
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x  +  y 


2^  1  dXi 

(A. 4) 


Now,  if  we  let 


and 


k  2tr 
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(A. 5) 
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then  the  resulting  integral  is 


(A. 6) 


wnere  the  constant  term  in  front  of  (A. 4)  has  been  dropped.  With  the  ad¬ 
ditional  definitions 


©  =  Vp"  V 
d©  =  dV 

the  definite  integral  of  (A. 6)  is  smitten 


(A. 7) 


oo 


O 


(A.e) 


(Ref.  13:380);  evaluation  of. the  integral  can  be  found  in  any  handbook  of 
mathematics.  Substituting  the  value  of  (3,  from  (A. 5),  into  (A. 8)  yields 
the  constant 


I  Ll 

V  2i§ 


With  the  integral  evaluated,  (A. 4)  becomes 


2g 


(A. 10) 


If  the  exponential  is  expanded  with  a  common  denominator  of  4Lg  ,  then 


'+9 


(A.10)  takes  the  final  form  of  the  kernel  when  the  substitution  of  2rr/X 
is  made  for  k,  i.e., 

exp{-  2^  -  !)(*»  +  yE)  -  2xy]}  (A.ll) 
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APPENDIX  B 


Before  equation  (2.3l)  takes  the  form  of  equation  (l.3)>  substitu¬ 
tion  of  some  variables  must  be  accomplished. 

-  ^(M  -  rr)y2  -iTfF[g(x2  +  y2)  -  2xy] 
g(y)  e  e  dy 

(B.i) 

Given 


_  _  ±\ys 

Vg(x)  e  2  ‘  =  -M 


"j 


(Ref.  11:15?)  and  solving  for  g,  we  get 


(3.2) 

(3.3) 

(B.4) 

(3.5) 

(B.6) 
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APPENDIX  C 


The  coefficients  of  the  polynomial  are  not  readily  determined  from 
the  form  listed  below;  however,  with  the  following  definitions  it  is  the 
most  concise  form  available. 


F  =  E  vN_n  F  (a) 
a  nv  ' 

Fp  =  £  VN'n  F  (g) 

P  n 

G  =  Z  vK_n  G  (a) 
a  n 

Gft  =  I  vN_n  G  (g) 
P  n 


(C.l) 


Here  it  is  understood  that  the  summation  is  from  n  =  1  to  n  =  N  . 
The  polynomial  is  now  written  as 


v 


2N+1 


V2”  -  ♦  Ga)  *  -  fh+1  ♦  Ga  -  C.J+1) 

+  '’<FeGa  -  F<xV  +  <FaG6  '  FgGa> 

+  Ft.+l(Ga  *  GS>  -  S+l<Fa  -  V  '  0 


(C.2) 


Also,  the  Fj.+^  and  G^.+^  3X6  the  series  functions  with  the  x  depen¬ 

dence  neglected. 

The  constants  a  and  b  are  then  calculated,  once  the  eigenval- 
n  n 

ues  have  been  determined,  by  specifying  a  particular  mode, 
equations  of  (^.21)  are  combined  to  yield 


The  first  two 
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APPENDIX  D 


The  plots  presented  here  are  a  few  of  the  higher  loss  modes,  and 
are  included  for  additional  verification  of  the  computer  program  by  com¬ 
parison  with  the  graphs  found  in  Reference  5*  All  but  the  last  two  plots 
are  a  result  of  the  second  order  approximation.  Figures  D-9  and  D-10 
are  seen  to  be  the  intensity  and  phase  across  the  feedback  mirror  only, 
and  axe  identical  in  structure  to  Figures  6-1  and  6-2. 
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Fig.  D-10.  Phase  plot  across  the  feedback  mirror  only  using 
the  first  order  approximation. 
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APPENDIX  E 


A  listing  of  program  BARC2  employing  the  analysis  and  derivations  of 
the  preceding  chapters  is  given  here.  Also  included  is  a  list  cf  input 
variables  required  for  program  operation.  Unless  otherwise  stated,  the 
input  variables  follow  normal  Fortran  convention  for  being  either  real  or 
integer  values. 


3ARC2  Inputs  (variables  are  listed  in  order  required)  - 

MAG:  Cavity  magnification  (real) 

NEQ:  Equivalent  Fresnel  number  (real) 

DELTA:  Mirror  misalignment  (offset  in  fraction  of  mirror  radius) 

NBIG:  Desired  number  of  terms  in  field  series 

MTEST1 :  Input  0  to  list  eigenvalues. 

MTEST2:  Input  0  to  continue  with  other  options,  1  to  do  new  cav¬ 
ity  (MAG,NE)5, DELTA) ,  or  2  to  exit. 

MODE:  Desired  mode  number  for  phase  and  intensity  calculations 

( 1  to  2>U3IG  +  l).  The  eigenvalues,  A.,  are  listed  accord¬ 
ing  to  loss,  i.e.,  mode  #1  is  the  lowest  loss  mode. 

MTE3T3:  Input  0  to  calculate  intensity  and  phase,  1  to  continue 
with  other  options,  or  2  to  exit. 

MTESTd;  Input  0  to  calculate  intensity  over  expanded  range.  If 
MTESTb  equals  0,  control  is  to  subroutine  ALLINT  with 
variables  II'CX  and  MTEST5  skipped. 

INCX:  Increment  value  of  x  for  phase  and  intensity  calculation 

MTEET5:  Input  0  to  list  field,  phase,  and  intensity  across  the 
feedback  mirror. 

MTEBT6:  Input  0  to  plot  constants  and  b  versus  NBIG,  1 

to  return,  or  2  to  exit  (return  is  LonMTE'T2).  If 
MTE3T6  equals  0,  the  plots  are  generated  and  control  is 
automatically  returned  to  MTE3T2. 

Mote  -  INCX  is  the  integer  number  of  points  between  consecutive  whole  x 
values.  For  example,  INCX  =  100  specifies  100  points  between  x  = 
0  and  x  =  1  for  the  phase  and  intensity  calculations. 
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Subroutine  ALLINT  Inputs  - 


■  XMIN:  Minimum  x  value  over  which  intensity  and  phase  are  cal¬ 
culated. 

XMAX:  Maximum  x  value  over  which  intensity  and  phase  are  cal¬ 
culated. 

IiJCX:  Increment  value  of  x 

NTEST1:  Input  0  to  list  field,  phase,  and  intensity. 

N  TEST  2:  Input  0  to  plot  intensity  and  phase. 

Program  operation  is  returned  to  MTEST6  in  BARC2. 


Note  -  For  the  best  results  let: 

and 


XMIN  =  -MAC  +  DELTA* MAG  -  0.5 
XMAX  =  +MAG  +  DELTA*MAG  +0.5 


-  Since  ZCPOLY  (the  root  finding  routine)  limits  the  degree  of  the 
polynomial  to  49,  the  maximum  value  for  NBIG  is  24  . 


100=  F'ROGR'AM  PARC? (DAT  A,  INCH  I  •  Oil  I  PUT  ,  TAPL  8=0inPIJT  .  T  AF'E  4  =  DA  I  A  ) 

110=  COMMON  SUTREX<70Q> »  PHASE ( 700 > . X INTFN < 700 ) * PHASE? <  700 ) 

120=  REAL  NE0»MAG.M0UPN<51 ) ,MSUPN<51 ) 

130=  COMPLEX  E  YE  »  COEF'  <  51  )  » LAMBDA  <  SI  )  *  CONS  I  A  <  SI  >  •fill  0(700) 

140=  COMPLEX  CL  <51  )  .CDUM.AN1  .AM2.X1  »  X?  »  Y  t ■  Y?  »  7.1 . 7. 2 

150=  COMPLEX  ROOT  »C0NSTD<51  )  ,PC0NA‘<51  )  .PCl!NB<51  )  .F  NX.CiNX 

160=  COMPLEX  r-  PEI  A  ( SI  )  *  F:  ALPHA  <51  >  *  GRETA  (51  )  »  GAI  PIIA<51> 

170=  DIMENSION  Fv'INIiF'X  <  51  )  .LABEL <25 > 

100=  HAT  A  LABE'L/25  ( 1  OH  )/ 

190  C*t****)**t ****************************  **************** ******* 

?oo=c 

:>10--C  THIS  PROGRAM  COMPUTES  RLSCINATOR  MODE  El OENUAl. LIES  AND 
220*C  SURGE DUE NT L Y  EVALUATES  INTENSITY  VALIIFS  EUR  POINTS 
230=C  ACROSS  THE  OUTPUT  PLANE  OF  A  STRIP  LASER  RESONATOR. 

24<)=C  THE  PROGRAM  DEALS  ONLY  WITH  A  DARE  CAVITY. 

250=  C  OUTPUT  CONSISTS  UP  AN  EIGENVALUE  LIST.  WITH  PHASE 
260=0  ANI)  MAGNITUDE.  FIELD  VALUES  FOR  A  SELECTED  MODE 
270  =  C  EITHER  ON  OR  OF  F  THE  MIRROR.  PLOTS  Or  r  3  FLU  SE  RIES 

2S0=C  FUNCTIONS  C)R  WEIGFTTINfi  CONSTANTS.  AND  PLOTS  OF  INTENSITY 

290  =  0  ACROSS  THE  OUTPUT  PLANE  WITH  EITHER  LIMITED  UR  EX  TENTED 

300=C  RANGE . 

31 0=C  COMPILED  CODE  NEEDED  AROUND  110000  OCTAL  TO  LOAD. 

320=C 

330  =  C  INPUT  (TUANT  I  TIES  ARE  AS  FOLLOWS! 

340-C  DELTA  =  MIRROR  MISALIGNMENT  (OFFSET  IN  FRACTION  OE  MIRROR  RADIUS) 

350=  C  MAC.  =  CAVITY  MAGNIFICATION 

360=C  NEQ  =  EELUI VALENT  FRESNEL  NUMBER 

370=C  NBIG  =  DESIRED  #  TERMS  IN  FIELD  SERIES 

3S0  =  C 

390=C  NOTE!  FVMAG  DF NOTES  EIGENVALUE  MAGNITUDE.  AND  F VPH  DENOTES 
4 0 0  =  C  EIGENVALUE  PHASE. 

410  =  1  THIS  PROGRAM  ALSO  REQUIRES  IMSt  ROUT  INI  ZOPOI.Y  AND  PLOTTING 
420=0  ROUTINE  CCPL0T56X.  FINAL  COPY.  20  SIFT  19B1  .  R  W  DERD1NF. 

430=C 

'I40--  C******  *******************************************  ************ 

450=  PI=2.*ASIN< 1 .0) 

460=  EYF  =CMf'LX  (  0  -  » 1  .  ) 

470  -994  FORMAT <  F  6 . 2  «  4X ) 

480=990  FORMAT  <  1X.F5.4.4X) 

490=  LAPEL  <  1  )  =  1  OFT  NEQ= 

500=  L.ADEL  <  3 )  =  10H  MAG  = 

5.10=  LABEL.  (13)  =1  OH  DELTA* 

520=777  WRITE < 0.999s 

530=999  FORMAT (1H1. IX. *  INPUT  MAG.  NEG .  AND  OFFSET!  *./> 

540=  READ  *»MAG» NEQ. DELTA 

550=  WR I I E <  8 • 00 ) MAG • NEQ »  DEL  T  A 

560=  WRI TE <  4 ) MAG . NEQ . DELTA 

570=  L ABEL. < 5 >  =  10H  MODE  RIG 

580=  LABEL  (  6  )  =  1  OF  IF  N  VALUE  ! 

59  ()~C** ***************** *****************  ************************* 

600=C 

6 1 0  =  C  MSUPN< I )=MAR**< I- 1 > 

620=C  MSUBN ( I ) = 1+ 1 /MAG# #2  +  ...  F 1/MAG** < 2* 1-2 > 

630=C 

A40=C****** *****************************************  ************** 

650=  MSUBN < 1  )  =  1 .0 

660=  MSLH’N  <1  )  =  1 . 0 

670=  DO  10  1=2.51 

6H0=  MSI  IF 'N  (  I  )  =  MAG*MSLIF'N<  I  —  1 ) 

690=  MSUBN <  I  )  =MSUBN <1-1)41  ./MSIJPNO  )**2 

700=10  CONTINUE 

710=  I  =2*NEQ*F' I  *MAG**2/  <  MAG**?  1.5 

72"  =  RNDIG=AL()G  <  ?50*NFN  )  /Al.PG  <  MAS ) 

730=  IF (RNBIG.LF  .50. >  GO  TO  11 
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740=  WRITE <8>998> 

750=  GO  ro  777 

7*0=11  URI TE  <  8 » 99A  )  RNF<  IG 

770=99*  FORMATE  IX.  #CALCUl.ATEli  NBIG  =  *  .  G1 4 . 7  •  *1  NflJT  INTEGER  CHOICE.**  •/) 
780=  READ  *  » N  F<  I G 

790=  WRI TL(8.979)NBIG 

800=  UR  I  FE  <  8 . 993 ) 

810-993  FORMAT < IX  .  *INF'UT  ZERO  TO  LIST  EIGENVALUES  :*./) 

820=  RE  All  *  »  MTEST1 

830=  WRITE  <  8  »  979  >MT  ESI  1 

840=C* ************************************** ********************** 

850=C 

8A0=C  COMPUTE  COEFFICIENTS  OF  THE  POLYNOMIAL 

870-f  F-(Z)=C0EF<1  )*Z**NMEG  -1  COEF  <  2 )  *Z**  <  NIiEG-1 )  +  ...  + 

GHO=C  COEF  <  NIiEG )  *Z  +  COEiF  (NLiEG+1  ) 

090=C 

900=C**************** ********************************************* 

910=  ALE HA=~1.+ DELTA 

920=  HETA=1 . +DELTA 

930=  coe:f<  1>=CMPLX<  1  .  »0.  ) 

940=  NBEG=2*NBIG-F  1 

950=  NC0EF=NDEG+1 

9A0  =  M=NP I G+ 1 

9/0=  HO  15  1=1 .M 

980=  AN1=CS0RT(4. *EYF*RI*T/MSUDN < I  )  ) 

990=  AN2=- T*EYE/MSUBN< I > 

1  000=  AN3=BETA*  <1.-1.  /MSUF'N  < I +1 >  > 

1010=  AN4=ALF'HA*  <1  .  - 1  .  /MSUF’N  <  1+  1  )  > 

1020=  AN5=BET  A-ALF'HA/MSUE'N  <  I  +  1 ) 

1030=  ANA=ALF'HA-BET  A/MSUF'N  <  HI) 

1040=  FBETA< I ) =- < CEXF < AN2*AN3**2 > /AN3 ) /AN1 

1050=  FALF-HAd  >  =- <  CEXF  <  AN2*AN5**2 ) /AN5 ) /AN1 

1060=  GF<£  I  A  <  1  )  =  <  CEXF-  <  AN2*ANA**2  >  /ANA  )  /AN1 

1 070=  GALF'HA  <  1  )  =  <  CEXF' <  AN2.*AN4**2  >  /AN4  )  / AN1 

1080=15  CONTINUE 

1  090=  FAL  F'HA  <  M  )  =CE‘XF'  <  AN2*BETA**2  > /BETA  /  <  -AN1  > 

1100=  FBETA  <  M  )  =  F ALPHA  <  M  ) 

1110=  GBETA  <  M  )  =CEXF'  <  AN2*Al  F‘HA**2  )  /ALF'HA/AN  1 

1120=  GALPHA(M>=GDETA(M> 

1  130=  COEF  <  2 )  =-  <  FBFZTA  (  1  >  +GALFHA  <  1 )  +  1 .  ) 

1140=  L1  =  NBIG  F2 

1150=  NA= 1 

1160=  NB=1 

1 170=  HO  21  1  =  3 » LI 

1180=  X2=CMF'LX  <  0  .  *  0  .  ) 

1190=  Y2=.X2 

1200=  LiU  18  JA=  1  ?  NA 

1210=  KA=NA- JA+ 1 

1  220=  X 1  =FDETA  <  JA  >  *GALF'HA  <  h'A  )  -FALF'HA  <  JA  )  *GBETA  <  K  A  > 

1230=  X2=X 1+X2 

1240=18  CONTINUE 

1250=  IF< I , EQ . 3 )  GO  TO  20 

1  260=  Ilfl  19  JB=1.NB 

1270=  KH=NB= JBT 1 

1280=  Y1=FAL  F'HA  <JB)*GBETA(KB)-FBETA<JB) *G ALPHA  <  KB ) 

1290=  T2=Y1+Y2 

1300=19  CONTINUE 

1310=  NB--NB+1 

1320=20  Z1  =F"BETA<  I -2) -FBETA <  1-1  ) 

1  330=  Z2=G AEF'HA  <  I  2  )  -PAl  F'HA  <  I  - 1  ) 

1340=  COE  F  < I ) =X2  F Y2+Z1 +Z2 

1 350=  NA=NA+1 

1 3A0=7 1  CONI  I NIJF 

1370=  L2  =  NHIG  E  3 
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1300=  NA=NBIG-1 

1 390=  NIi=NBIG-l 

i4oo-  no  :»o  i=l2.ncoee 

1410=  X2=CMPLX  <  0  .  *  O  .  ) 

1420  =  Y2=X2 

1430=  lTd  .EH.NPUEF)  GO  10  23 

1440=  no  22  JA=1 .  NA 

1400=  KA=NB)G1  JA-NA 

1460=  X  1=EI<ET  A(  M- JA )  &GALPHA  ( KA ) -FAL  I  HA<  M- JA ) #GBET  A ( KA ) 

1470=  X2=X 1 4X2 

1400=22  CONT INUE 

1470=  ird.E0.L2)  GO  TO  7T. 

1500  =  23  DO  24  JF<  =  1.NB 

1510=  KB=NBIGf JB-NB 

1  520=  Y 1  =F  ALPHA  <  h  -  JB  )  *GBETA  <  KF  >-FBET  A  <  M- JB  )  *  GAL  PH  A  ( KB ) 

1530=  Y2=Y 1  FY2 

1540=24  CONTI  NIJF. 

1550=  NF=NF-1 

1560=  GO  TO  27 

1570=25  HO  26  .JE<=1.NPIG 

1 580=  KF=NBIG  JBH 

1  590=  Y  1  =  1  ALPHA  <  JB )  *GBETA  <  KF  )  -EBET  A  (  JF  )  *GALPHA  (  KTi ) 

1600=  Y2=Y 1 +Y2 

1610=26  CONTINUE 

1620=27  Z 1 =EBET A  ( M  )  *<  G ALPHA  < I  -M-l > -  DUET  A  <  I ~M  1 ) ) 

1630=  Z2=GALPHA  <  M  )  #  <  F  ALPHA  <  I  -M-  1  )  -F  BF-T  A  (  I  -M  1  >  > 

1640=  COEF  < I ) =X2+ Y2 f Z1 -Z2 

1650=  NA=NA- 1 

.1  660  =  2EI  CONTINUE 

16  70=0********  ************#*******.  I-********************?*  ********** 

1 600=C 

1 6 9 0  =  c  COMPUTE  K-OOTG  OF  POLYNOMIAL  WITH  IMSl.  ROUTINE  ZCPtlLY.  TO 

1 700 =C  OPT  AIN  THE  EIGENVALUES.  AND  THEN  ORDER  EIGENVALUES  FT 

1710=r;  SIZE. 

1 720=C 

1. 730=C******************  ***********************  *********  ****  ******* 
1740=  CALL  ZCPOLY < COEF rNDEGr LAMBDA. IER> 

1750=  IF (MTEPTt . EO.O)  WRITE <8. 8V) 

1760=  1=1 

1770=  DO  70  T  1-2  r  NTUIEi 

1780=  !3 1  ZE=  REAL  <  LAMBDA <  I  >  )  **2TAIMAIT  ( LAMBDA  ( I )  >**2 

1790=  K  =  I 

1800=  DO  75  J=:t  t.NDEG 

1810=  SIZE) -  REAL (LAMBDA! J) ) **2+AI MAG ( LAMBDA ( J > )**2 

1820=  IFCSIZEl  .LT.SI7-E)  GO  TO  75 

1830=  k=J 

1040=  SI 7E  =  ST  ZE 1 

1.850  =  75  CONTINUE 

1060=  CDUM=l AMBDA < I ) 

1 870=  LAMBDA < I )= L AMBDA < L ) 

1880=  LAMBDA  ( K  )  =  CDIIM 

1890=  CL ( I ) “LAMBDA ( 1 ) 

1900=  EVPH=AT AN2<AJMAG<CL<  I ) ) »RFAL(CL< I > ) ) *180. /PI 

1  910=  SMA=REAl  <  CL  <  I  )  >  **'.M  AI  MAG  <  CL  <  I  >  >**2 

1920=  SMAG  =  SOR T  <  SMA ) 

1930-  IF (MTESI 1 .EU.O)  WRITE ( 8. 333 ) I » l AMBDA < I ) » SMAG . FVPH 

1940=333  F  ORMAT < 1 X , I  9 . 1 X . 4 < G1 4 . 7 . IX >  »  7  ) 

1950=  1=11 

1960=70  CONTINUF 

1970=  EVPH*A  T  AN2  <  AII1AG  <  LAMBDA  (  NIlEG )  >  .REAL  ( l.AMBDA<  NDE  G  >  )  >*180. 7  PI 

1980=  SMA=F<EAL(  LAMBDA  <NHEG)  )  **2  +  A IMAU  (  LAMBDA  <  NTiEG  )  >**2 

1990=  SMAG=SHRT (SMA) 

1000=  IF  <  M  T  LSI  1 . L 0.0)  WR I TE < 8 . 333 >NDfc 6 . 1 AMBDA  <  NIC fi ) .SMAG.LVr H 


r-j  rj  rj  ro  ,j  rj  rj  u  rj  rj  ,j  rj  rj  rJ  rj  rj  ,j  rj  »o  rj  rj  ,  j  rj  (J  fj  r:  r-j  rj  fj  rj  tj  ■  j  f  j  ro  f_'  rj  rj  rj  r_'  f j  f j  rj  ,-j  rj  r J 


=C**  **************************************************  ************ 

=C 

■=c  NOW  calculate:  the  constants  for  the.  function  sum  tor  a  par- 

=C.  T  I CUt  AFC  MODE.  THEN  LOOP  TO  CALCULATE  T  MF.  TIKLH  A  T  A  SEIECUH 
-C  NUMBER  OF  POINTS. 


i'************************************************ ****** ********* 


TO  HO  NEW  CAVITY, 


WRI T E  < 8 ,997 ) 

FORMAT  ( IX,  *  INPUT  0  TO  CONTINUE,  1  TO  HO  NEW  CAVITY,  OR  2  TO  E 
I  ,/) 

REAR  *»MTEST2 
WRITE  <  8 , 97V  )  MTEGT2 
1 F ( HTEST2- 1 )  100,777.088 
WRI  TE  <  0 , 990 ) 

T  ORMAT  (  1  X  ,  *  I NPU  T  HES1REH  HO  HE  NUMEiF.R  L  *  ,  /  ) 

REAR  *.MOHE 

WRI  I E  ( 8  •  97V  )  MORE 

LABEL  ( 1 5 )  -  1  OLE  MODF  * 

ENCODE  < 10 « 9B7, LABEL < 16) ) MODE 
FORMAT (4X» T  2 , A  X ) 

ROOT  =LAMRHA <  MODE: ) 

WRIT E< 4) MODE. ROOT 
X2=CMPL  X  <  0  .  ,0.  ) 

Y2  =  X2 
Z2  =  X2 

no  40  i  =  i,Nnic, 

RTNHEX ( I )-I 

AN  1  "ROOT ** ( NRIG-I  > 

X1=AN1*TALPHA(I )*<R00T-1. >/FBETA(M) 

X2=X 1 +X2 
Y1=AN1*GALPHA( I) 

Y2  =  Y1  EY2 

.  Z1=AN1*T  ALPHA ( I ' *GDETA ( M >/F  BETA ( M > 

Z2-Z1 f Z2 
CONTINUF 
HO  41  I-1,NHIG 
AN1=R00T**(NBIG- I ) 

CONSTB  <  I  >=AN1*(X2i 1  .  )  /  ( ROOT**ND IG-Y2F7.2 ) 

CONST  A( I ) ~ ( AN  1* (  ROUT-1 .  > "CONST B( 1 >*GHET A(M) ) /FRET A < M > 

CONTINUE 
WRITE:  (8, 982) 

FORMAT  <  IX. *INPIJT  0  TO  CALC  INTENSITY,  1  TO  CONTINUE,  OR  2  10 
1  *./) 

READ  *, MI  ESI  3 
WRITE ( 8 , 979 ) MTE  S I  3 
ENCODED 10, "94.1  ABEL (2) )NEO 
ENCODE < 10, 994, LABEL (4 ) ) MAG 
ENCODE (10«  990 , LABEL (14)) DELTA 
ENCODE < 1 0,980. LAHEL<7> ) REAL (ROOT ) 

ENCODE  (  10, 980 . t. ABEL  ( 8  )  >  AIMAG ( ROO  I  > 

FORMA  T ( 010 . 4 ) 

TF <  MT  ESI  3- 1 )  101,102,880 
t  ABEL  (  9  )  ~1  OHM  J.RROR 
t  ABEL.  (10  )  =  1()HTT.  ANE 
LABEL  < 11 >  =  10HSCALEH 
LABEL( 12)=10H  INTENSITY 
LABEL ( 1 7 )  =  1  OH  PHASE 
LADE l  ( 18) "  I  OH  (DEG) 

WR I TL ( 8 , 222 ) 

FORMAT ( IX, *INPOT  0  TO  CALC  INTENSTIY  OVER  EXPANDED  RANGEL*./) 

READ  * , MTEST  4 

WRITE(8,97?)MTF.ST  4 

IF  <  MTEST  4 . NL . 0 )  GO  TO  104 


TO  E  XJ  T  :  * 


INTENSITY, 


TO  CONTI NOE, 


2  10  EXIT 


. 


2640  = 
2650  = 
2660  = 
2670= 
2680  = 
2690  = 
2700= 
2710  = 
2720  = 
2730= 
2  740= 
2750= 
2760  = 
2770  = 
2780= 
2"790= 
2800= 
2810  = 
2820= 
2830  = 
2040  = 
2850= 
2860  = 
2870  = 
2880- 
2890  = 
2900  = 
2910  = 
2920  = 
2930  = 
2940  = 
2950= 
2960  = 
2970= 
2980 
2990  = 
5000 
3010 
3020 
3030 
3040 
3050 
3060  = 
3070  = 
3080= 
3090= 
3100  = 
3110= 
3120= 
3130= 
3140  = 
3150  = 
3160  = 
3170 
3180  = 
3190  = 
3200  = 
3210  = 
3220= 
3230 
3240 
3250 
3260 
3270= 


CAL  L  ALL  1  NT  (  MAG  .  MSI1BN  .  MSI  IV  II  .  CONST  A  .  CONS  I  8  .  T  .  NLUG  •  LOU  I  , 

1  ALPHA. BF TA.t ABEL . HCLTA » MODE  ) 

GO  TO  102 
104  WRITE  <  8 .992 ) 

992  FORMA  MIX.  *  INPUT  INCREMENT  Of  X  FOR  PLOT.#*/) 

REAti  *  ,  I N  C  X 
WRIT! (8,979>INCX 
WRITE  (8,950) 

=950  FORMAT (IX .* INPUT  0  TO  LIST  FIELD.  PHASE.  ANH  INTENSITY  :  *./) 

=  RE  AH  *  »  MT  EST'.i 

=  WRI I E  ( ft .  979  TMTEST  5 

=  I F ( MTEST5 . NE . 0 )  GO  TO  556 

WRITE (8. 554) 

554  FORMAT  <  17X  »  *FTELD*  »  20X »« INTENSITY*  .  1 IX  ,  *X*»  11X»*PHASE  (  DEG  )*»/') 

556  NHATA=0 

X - ALPHA 
IiRIGHT  =  0. 

31  NHAlA=NHATAt 1 
ST  OREX ( NDAT  A ) “X 
X 1 =CMPLX ( 0 • .0. ) 

X2=X  1 

DU  32  J  ==1  »ND1G 

AN  1 =CSQR  t <4.*EYE*PI*T/MSUBN< I ) ) 

AN2=-  T  *1  YE/MStlfiN  (  1  ) 

AN3=ALE'HA-  X/MSUf  'N  (  I  +  1  ) 

AN4-BF  T A  -X/MSUPN  (HI  ) 

ENX=  CONST  A ( I >*CCXP<AN2*AN4**2)/AN1 /AN4 
GNX==C0NSTB<I)*CEXP<AN2*AN3**2)/AN1/AN3 

XI- -FNXFX1 

32  X2=GNX+X2 

FIELD (NDATA)“1+X1+X2 

XJN  T  EN  <  NIlAT  A )  =  REAL  (f  IEl  II  (  NDATA  )  )  **24  A1  MAG  ( F  I  El  DUIDATA  >  )**2 
F:,HASE:  (  NI.i  A I A ) «  AT  AN2  ( A I  MAG  ( F I  ELII  (NIiAT  A ) )  » REAL.  ( F  IELD(NDATA)  >  ) 

PHASE  <  NliAT  A  )  = PHASE  (  NBA  TA  )*1 80  .  /PI 
AT  A  1MAG<  I  lELLKNDAIA)  ) 

A2  =  -REAl.<  FIELD  (NDATA)  ) 

PHASE 2  (  NDATA  )==ATAN2  ( A1 »  A2  ) 

PHASE2 ( NDAT A ) =PHASE2  <  NDAT  A ) *180 . /F  I 

IF  (  X.TNTCN  ( NDAT  A  >  .GI  .  BRIGHT  )  BRIGHT =X INI EN < NDAT A > 

IF (MTESTS.NE.O)  GO  III  557 

WRITE  <  8 . 555 )F IELH ( NDATA  ).X1N  ILN(MMATA)  *  STORLX  (NDA >  A )  •  PHASE  ( NHAT  A  ) 

555  FORMAT <  5X . 2G1 4 . 7 . 5X . G 1 4 . / . 5X . fil 4 . 7 . 5X . F7 . 2 ) 

557  X  =  XF  1  .  /IN< -X 
TF(X.LE.DETA)  GO  TO  31 
N1=NI'ATA+1 

N2=0 

PHASE < N1 >=0. 

HO  37  I-l.NHATA 
I  MET  A=aHS  PHASE  <  M  1  )  -PHASE  <  I  >  > 

IF  (  THETA,  (ir  .300.  >  N7=N2il 

37  X I N  TEN ( I )  =  X I N T EN ( I ) /HR 1 GHT 

CAL  t  HGRAF'H  ( ST  OPEX  .  X 1 N I  EN .  NDAT  A .  I  ABEL.  .1.0.0) 

IE<N2.GT.1>  GO  TO  38 

CAt  L  HHRAPH  (  S  TURF.X.J  PHASE.  NHAT  A.  LABEL  ,3.0.0) 

GO  TO  102 

38  CALL  FIGRAPH  <  ST  OREX  ,  PHASE.?. ,  NDAT  A  » l.  ABF  L  .  3 . 0 , 0  > 

102  WRIT  E  <  G » 991 ) 

991  FORMAT!  IX. *IYPE  0  TO  PLOI  CONSIS  OS  N»  1  It)  Rl  TURN,  OR  2  10  EXIT! 
1  *./> 

READ  *.MIFST6 

WRITF(Ll,V79)MI E  ST 6 

IF ( MT  ESI  6  1)  103-45.888 

103  I  ABEL (9)-10HCUNSIANT  * 

=  I  ABEL ( 1 0 )  =  1  OH 
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3200=  LAPEL <11  )  =  1  OH  MOmCUNST 

3290=  L ADFL < 1 2 ) =1 OHANT-  A>**2 

3300=  DU  42  I'liNUG 

3310=  F'CONA  <  1  )  =RLAL  (CONST A <  I  )  )  **2  t  AIMAG  (  CONS  1 A  <  I  )  )  **  2 

3320=42  PCONIK  I  )  =RE Al.  (CUNSTpF I  )  ) * *2  FA I MAG < CONST D (  1 )  >**2 

3330=  CALL  HGRAF'H  ( R INHEX  •  PCONA . NDTG  »  LAPEL  »  1  »-l  ,  1  1  ) 

3340=  LAPEL ( 1 2 )= 1 OHANT  10**2 

3350=  CALL  HURAPH  <  R  1  NDEX  *  F'CUNP  •  NPI  (3  ,  L.ABFl  ,1,-1 ,11) 

3360=  WR] I fc ( 0* 904 ) MLIDE 

3370=984  FORMAT  (IX*  *CUMF'LET  I  D  PLOT  OF  CONSTANTS*  HOOF.  =*,12,7) 

3380=  GO  TO  45 

3390=970  FURMAT  <1X,*INP(JT  VALUE'S  ARE  5  *,2I5,/> 

3400=979  FORMAT < IX, *  INPUT  VALUE  IS  !  *,15,/) 

3410  =  998  FORMAT <1X,*REV1SL  PARAMETERS  SO  NDEG  <  24*) 

3420=88  FORMAT  <  1()X,*(1AG  =  *  »F6 . 3 , 5X  ,  *NEP  =  *  ,  I  6 . 3 »5X  ,  *HEl  I A  =  *, 

3430=  1  F6.5,/) 

3440-89  F  ORMA !(9X,*I*,2X, CLAMBDA <  REAL ) * , 3X , *L AMDDA < I MAG  >  * , 6X , 

3450=  1  *E."  VMALF*  ,  1  OX  >  *EVPH*  ,  /  > 

3460=888  CALI.  EXIT 
3470=  ENU 

3480=  SOPROOT  TNE  All  J  NT  <  MAG  ,  MSUBN,  MSUF'N  ,  CONST  A  ,  CONS  I  P  »  T  ,NPIU,RUOT  , 

3490=  1  ALF'HA  ,  T'E  T  A  ,  L  APEL ,  PEL  T  A  ,  MOPE  ) 

3500=  PI MENS TON  L  APEL ( 25 ) , NCOUNT  < 1 0  > , J  <  25 ) 

3510=  COMMON  S10REX<7<)0>  ,  F'FTASE  <  700  )  ,X INTENT  700)  , PHASE? v  700) 

3520=  RFAl.  MSIIPN  <  51  >  ,  MSIJF'N  <  51  )  ,  INARG1  •  I NARG?  •  I  NARG3  •  INARG4.  FNAR05 

3530=  RF AL  INARG6 , LNARG7 , INARG8 , MAG  *M ] NO 

3540=  COMPLEX  APART  1  ,  APART2  ,  PPART  1  ,  KPAR  T  2  ,  ALL  FUN  , UIHS  T  A <  5 1  )  .ROOT  ,l'.Yt 

3550=  COMPLEX  AFUN ,  PF  UN  •  SPNTUO  ,  SENT  111  ,  1)0*  ,  OUT  CON ,  FRESL  •  CONST  P  (  51  >  • 

3560=  COMPLEX  E YEFAC ,EYEF 1 ,EYEF2 , SPNTU2 , VAR4 ,r IEl  IK  700 )  • 

3570=C**********»*** *********************************************** ********** 

3580=C 

3590=C  THIS  SUDKOUT INE  F  OLLOWS  PROGRAM  PARC  ANP  COMPUTER  REAM  INTENT! IT  1 T f  S 
3600=C  TN  THE  OUTPUT  PLANE.  *  INTERMEDIATE  POINTS  LOR  EVALIIAUON  ARE  INPUT 
36 1 0=C  WHILE  ALL  OTHF-R  REUUIREP  QUANTITIES  ARE  CARRIED  I H ROUGH  IN  THE 

3620=0  ARPUMLNI  LIST  AS  FOLLOWS.4 
3630=C 

3640=C  MAG=  CAVITY  MAGNIFICATION 

36S0*C  MSUPN-  ARRAY  FOR  PARTIAL  SUMS!  OF  INVERSE  POWERS  OF  MAG 
3660=C  MSIIPN*  ARRAY  FOR  MAG  TO  SOME  POWER 

3670=C  CONST-  ARRAY  OP  CONSTANTS  IN  THE  ASYMPTOTIC  SERIES 

36B0=C  T=  QUANTITY  DEFINED  IN  PARC  PER  HORWIT  7 

3690 -C  NDIG*  I  TERMS  IN  THE  SERIES 

3700--C  ROOT*  MODE  EIGENVALUE 

37 1 0=C  LAPEL*  PLOT  LABELING  ARRAY 

3720=C 

3 730=C********* ******************************************************** ****** 
3740=  NDAT  A  =  1 

3750=  PRIGHT=0. 

3760=  SX 1=0.0 

3770=  SX2=0 . 0 

3780=  SY 1=0.0 

3790=  5XY=0.0 

3800=  DO  5  1=1,10 

3010=5  NCOON I ( T ) =0 

3U20=  XI =<-l ,+DFl TA'*hAG 

3830=  X?  =  <  1  .  TPF  1.  IA)*MAG 

3840=  F 1=2. *AS IN ( 1 . ) 

3050=  EYE=CMPl.X  ( 0  •  ,  1  .  ) 

3860=  F  YEF  AC=  < 1 . -F  YE ) /2 . 

3870=  WRITE (8, 900) 

3800=900  FORMAT  (1H1  .CENTERING  EXTENFiED  RANFiE  INTENSITY  SUBROUTINE.*,/) 

3890=  DO  10  t  =  1 .51 

3900=10  MSIJF'N  ( I )  =  MAO*MSIJPN  <  I  ) 

3910-  WRITE (8.901 ) 


3920=901  FORMAT  (  1  X  > » 1  Nt-’IJT  MIN  ANR  MAX  X  VALUES  ANR  «  POINTS  BETWEEN:  » 

3930=  Kt.AH  ♦.XM1N.XMAX.1NCX 

39-10=  WRITE  (R.9021XMIN.XMAX.  1NCX 

3950=902  t  LIRMAT  (  IX.*  INPUT  VALUES  ARE  :  *  *  F  ! . .  2  »  2  X  » F  5  .  •  ?X  *  1 5  »  /  ) 

3960=903  RIRMAT < IX.* INPUT  VALUE  IS  J  *,T5./> 

3970=  WRITE (8.904) 

3980=90-1  FORMAT  <  IX- *  INPUT  0  TO  LIST  F1LLR.  PHASE  •  ANli  INTFNSJTY  !  *./> 

3 990=  RFAIt  *.  NIEST1 

-1000=  WRITE  (  fl » 903  )  NTtSTl 

-1010=  If  (NTt  l.Tl  .NF:  .0)  00  TO  20 

4070=  WRITE (8. 90S) 

4030  =  900  FORMAT <  1  7X  .  *F I ELM* .  20X . *  I NT ENS1 T Y* , 1 IX. *X* . I IX . *PHASE  ( HE  0 ) * . 

4040=  1  11X.*PHASE2  <HEG)*./> 

4050=20  X--XMJN 

4060=50  ALL  F  UN  = (0 . • 0 .  ) 

4070=  HO  310  1=1. NR1G 

4080=  MFNV  =  1  ./MSLIPNC  I  ) 

4090=  VAR1  =2  .  *  <  1  .  +1  .  /MSLIPN  <  2#  I  > /MSURN (  1  )  ) 

4100=  VAR2=S0RT<4.*T/P1/VAR1 ) 

4110=  VAR3-MINV/MSUHNC 1 ) 

4120=  STAPHA=(X/MAGHtriA*VAR3)/< .5*VAR1  ) 

4130=  STAPHR=<X/MAG+ALPHA*VAR3>  /<  .5*VAR1  ) 

4140=  INARR1  =  <HETA-X/MAG )**:.*+<  RET A-r<ETA*MlNV>**2/MSUHN<  I  ) 

4150=  TNARG2=<RETA-X/MAG-RETA*(1 . - MINV ) *VAR3 > **2/ < .5*VAK1 ) 

4160=  A ARG 1  =  1 NARG2- T  NARG1 

4170=  APART 1=CEXP(E YE>T *AARG1 )* (-CONST A< I ) )/( RET  A- MET A*M1 NV ) 

1180=  INARG3-  < ALPHA-X/MAG)**2F < BET A-ALPHA*M1NV ) **2/MSURN < I ) 

4190=  1 NAR04= <  ALPHA-X /MAO- ( HET  A- ALPHA*M 1 NV  >  *VaR3 ) **  2/ ( .5*VAR1 ) 

4200=  AARG2=INARG4- 1NARG3 

4210=  APART2*CEXP < EYE*T*AARG2 ) *  < -CONSTA ( I ) ) / ( RET  A- ALPHA* M 1 NV ) 

4220=  J  NARGS=  <  HET  A-X/MAG )  **?♦  <  Al.PHA-PE  I  A*M1  NV  )  **2/MS(.IRN  ( I  1 

4230=  1NARG6= <  RE  T  A-X/MAG- ( ALPHA-RETA*M 1NV ) *VAR3  >  **7 / ( . 5*VAP 1  ) 

4240=  R ARG 1=1 N ARU6 - 1 NARG5 

4 250=  OPART  1  =CTXP <  BARG  1*PYE*T>*< CONSTB <!)>/<  ALPHA -IC IAJTM1NV) 

4260=  INARG7=  <ALPf  I  A-X/MAG  )**2T  (  ALPHA- Al.PHAfM.TNV )  *»2/MS(.)RN(  1  ) 

4270=  J NAMUR*  < ALPHA- X/MAG-ALPHA* < 1 . -MINV >*VAR3 > **2/< .5*VAR1  ) 

4280=  HARG2= I NARG8-T NARG7 

4290=  HPART2  =Ci:xP  (  RAKG7*EYE*T  )  *C0NS1  l»<I)/<  AL.PHA-AL  PHA*M1NV  ) 

4300=  (HIT  CON-SURT  <  MSLIPN  <  1  ) /PI/T /VAR1 ) /2 . /ROOT 

4310=  PRSA1  =VAR2*(HE  I A-X/MAG-RE 1  A*  ( 1 .  -MTNV)#VAR3i 

4  320=  F  RSA2=VAR2*  <  ALPHA-X /MAG-  <  RET  A- ALF'HA*M  t  NV )  *  VAR'3  ) 

4330=  FKSR1 =VAR2# ( RET A-  X/MAG - <  ALPHA-RET A*M 1 NV ) *VAR3 ) 

4340=  FKSR2=VAR2*< ALPHA-X /MAG -ALPHA* ( 1 . -MINV ) *VAR3 ) 

4350=  EYEF1=EYEFAL 

4360-  EYEF"2=EYt  F'AC 

4370=  I F < SI AF  FIA .01 .BETA)  EYEF2=~EYf FAC 

4380=  11  (STAPHA.GT .RETA)  EYEF  l-  EYEF AC 

4  390-  IF  (STAPHA.PE.Air-HA.  ANH.STAPHA.ee.  RET  A)  EYEF  2  - 1  YF  F  AC 

4400=  AMIN--APAR!  1  t(f  RF  SI  (F  RSA1  > -EYEF 1  ) ~ APART  2* (FRESH  (FRSA2)-F  YF'F2> 

4410=  SPNl<Il=CEXPl-EYE*(Pl/4.  +  T*< < ST APHA-X/MAG ) » *2 F ( RET  A- 

4420“  1  ST  AFHA*MINV)**2/MSI.IRN(  I )  )  )  )/<Rf  T  A- S  TAPHAtM  1  NV  1  *  <  -CONSTAT  I  ' 

4  430=  IF  (STAPHA.GF  .  ALPHA.  ANH.ST  APHA.LE.  RE  T  A)  AMIN=A1  (IN-FSPNT  HI  *GHRT  ( 

4410=  EYEF 1=EYEF AC 

4450=  EYEF  2=F  Yl  T  AC 

4460=  IF  (STAPHT  .GT.Rt.TA)  F.Y  E  F  l  =- E  YEF  AC. 

44/0  J  F  (SI  AF'FTR .01  .  F<F.  I A  )  F YEF 2= - EYEF AC 

44R0=  IF (  S  I  APFIFi .  GE  .  At  F’HA .  ANIl .  STAF’HH .  I  F  .  RT  1 A  )  EYEF 2  FYI'FAC 

4490=  RF  t)N=RFAK  I  1  *<F  RESl  T  F  RSR1  )  -  E  YEF  1  i - REAR  T  2*  ( F  Rt  51  <1  PGRI!)  I_  YEF'.?  ) 

4500=  SPNTFI2=CF:XP(-F  YE*<r  l/4.-TT#(  (STAPFIR-X/MAG)**2F  (ALPHA  5TAPHR* 

4510=  1  MINV)**2/MS(IRH<  I  )'))/(  ALPHA  S  I  APFIR*M  INV  )  *CONS  T  R(  1  ' 

4  520=  IF  (ST  ATT  IB  .  C»r  .  AL  F'HA  .  ANH  .  STAPFIR  .  l.fc.  .  RE  T  A  )  RFIIN=HF  UNI  SPNT  H7»S0H'T  ( 

4530-  ALL  FUN=(KI  ICON*  (  AFIIN  FRFUN  )  +  ALLFlIN 

4540=310  CUN1INUE 
4550=  EYEF  1  =EYFF  AC' 
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4560= 

4570= 

4500= 

4590= 

4600= 

4610  = 

4620= 

4630  = 

4640= 

4650  = 

4660= 

4670  = 

4  680  = 

4690= 

4700  = 

4710  = 

4720  = 

4730= 

4  740  = 

4750= 

4  760= 

4770= 

4780= 

4790  = 

4800= 

4810  = 

4820  = 

4830  = 

4840= 

4850  = 

4860= 

4870= 

4880= 

4890= 

4900= 

4910= 

4920= 

4930= 

4940= 

4950= 

4960  = 

4970= 

4»80  = 

4990  = 

5000= 

5010= 

5020  = 

5030  = 

5040  = 

5050= 

5060  = 

50  70  = 

5080= 

5090= 

5100 

5110  = 

5120- 

5130 

5140 

5150  = 

5160 

5170  = 

5180  = 

5190  = 


1  ^ 


FYEF2=EYEF  AC. 

11  (X/MAG.UT  .BETA,)  EYEF1=  eyefac 
IF  (  X/MAG  .  GT  .  BF.T  A  )  EYE  F  2=- L7EFAC. 

-  IF  <  X/MAG.GE  .  ALPHA.  ANU.  X/MAG.LE .  BE  I  A)  F  YEF2---F  YF.  1  A!. 
ERSU01=S0RT  (2.YT/P1  )  *  <  BET  A-X/MAG  ) 

FRSU02=SUR1  <  2  .  *  I  /F  I )  *  ( ALPHA-X/HAO  > 

VAR4=CSURI (EYE/2. >/R00T 

U0X=VAR4*<  <E  RESLtFRSUOl  )-EYEF  1  )-  (  FRESl.  (F  KSU02 )  -EYEF2 1  ) 

SI  NtU0=CEXF'(  -EYE  *PI/4  .  )  *VAR4*S(JRT  (  2  .  > 

I E  (  X/MAG  .  GE  .  At  PHA .  ANU .  X/MA  G  .  LE  .  BETA  >  U0X=U0X-1  SEN  1 110 
ALLFUN=ALLF  UN-E  UOX 
F  IEI.H  <  Nil  AT  A )  =-Al  LFUN 
STORF.X  <  NBA  I A  >  =X 

XINTEN(NUATA)=AIMAG<  ALL FUN  )#*2+E?EAL  (  ALL  F  UN  1**2 
F'HOSF.  (  NLiAT  A  )  =  AT  AN2  <  A I  HAG  <  AL  LF  UN  >  .REAL  (  Al  LF  UN  )  ) 

F-HASE  <  NBA  I  A  ) -PHASE  (  NMAT  A  >  *  I  80  .  /F'l 

RHASE2  <  NUAT  A  >  =ATAN2  < -Al MAG ( ALLFUN ) . -  REAL ( ALt  FUN  > > 

FHASE2  <  NBA  t  A  )  ==PHASE2  (  NUAT  A  >  #  1  80  .  /!'  I 

IF  <  XINTEN  (  NBA  t  A )  .  f!T  .BRIGHT)  BRIGHT =X1NT  EN  ( NUA I A  > 

IF  (  NTE.ST  1  •  NE  .  0  )  GO  TO  400 

WRITE  (8.906)  ALL.F  UN »  XINT  E  N (NUAT  A )  » ST OREX < NUAT A >  .CHASE:  <NHA1  A)  . 
1  PHASE  2 ( NBA  I  A ) 

906  F  ORMAT ( 5X • 2G1 4 . 7 . 5X , G1 4 . 7 . 5X , G 1 4 . 7 . 5X . F  7 . 2 . 1 5X . F7 . 2  > 

400  X=X+1 .  / INCX 

IF <X.GT .XMAX)  GO  TO  500 
NUAT  A==NI'A  I  A+  1 
GO  TF)  50 
500  N 1  '-NUAT  A+  1 

PHASE < N1 )~0. 

DO  510  I  ==  1  •  NUAT  A 
IF  <  STORE:  X  (  I  )  .LE  .XI  )  J(24)=I+1 
IF  < ST QRFX  ( 1  )  .LE  .X2)  .l<25)=I  +  l 
510  XI  NT  EN  < I )=XINTEN< I ) /BRIGHT 
K=-0 

.11  =.)(  24) 

12- J <  25 ) 

PCI  520  T  —  .11  » J2 

T  MET  A2~  ABS  < PHASES f  I  + 1 )  -PHASES <  I  )  > 

THETA=AHS(  PHASE  <1  +  1  )-PHASEO  >  ) 

I E  (  THETA2  .  GT  .  300 .  )  NCIOIINT  <  2 )  =NCGUN  I  (2  HI 
IE< THETA. LE. 300. )  GO  TO  520 
K=K  + 1 
J  ( K  >  =  T 

NCOUNT  (  1  )  =N0 CIUNT  <  1  )  + 1 
520  CUNTINUr 

UU  530  1=1.51 

530  MSUPN <  1 1 -MSUF’N <  I )  /MAG 

UFv'l'TE  C  8 » 907 ) 

907  FORMAT ( 1 X.» INPUT  0  TO  PLOT  INTENSITY  AND  PHASE  t  *•/> 

REAU  *»NTEST2 

WRITE <8. 903)  NTESF2 

IE  <  NTES I  2 . NE . 0 )  GO  TO  550 

CALL  HGRAPEKSTOREX. XINTEN. NtiAT A.LABEl  >1.0.0) 

IF  l  NCF1UNT  (  J  )  .  GT  .  1  )  GO  TO  535 

CAL  t.  HGKAPHt  STORE  X .  PHASE  .  NUAT  A .  LABEL  .3.0.0) 

GO  TO  550 
535  CONTINUE 

IF  <  NCUUN  T  <  2 ) . GT . 1  )  GO  TO  53  7 

CAL l  HPRAPH ( S 1 OREX . PHASE2 .NUA  T  A .LABEL .3.0.0) 

GO  TO  550 
5.3  7  k  =  1 

540  J 1  =  J ( K ) 

J2= J 1 K4 1 )-l 

=  if  <t\  .EO.NCOUNT  (  1 )  )  J2-J(25> 
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5200= 

5210- 

5220 

5230= 

5240  = 

5250= 

5260  = 

52  70  = 

5280= 

5290  = 

5300= 

5310= 

5320  = 

5330= 

5340= 

5350= 

5360= 

5370  = 

5380= 

5390= 

5400  = 

5410  = 

5420  = 

5430= 

5440= 

5450  = 

5460  = 

54  70= 

5480- 

54‘?0= 

5500  = 

5510  = 

5570  = 

5530  = 

5540= 

5550= 

5560= 

5570= 

5580= 

5590  = 

5600= 

5610= 

5620  = 

5630= 

5640- 

5650= 

5660  = 

5670  = 

5680= 

5690  = 

5700= 

571  0= 

5720  = 

5730  = 

5740  = 

5750  = 

5760= 

5770- 

5/80  = 

5790  = 

5800  = 

5810= 

5820= 

5830  = 


HO  545  1=J1,J2  ’ 

IF  <  PHASE  (  J1  )  .GT.O.O)  PHASE  (HI  >=PHASE<  1  +  1  >  +  36.0. 

If  (PHASE!  J1  >  .LT  .0.0)  PHASE  (  1  +  1  )=PHASl.  (HI  >-360. 

545  CONTINUE 
K=M2 

IF ( K • LE  .  NCCJUNT  ( 1  >  )  GC>  TO  540 
CALL  HGKAPH ( STOREX ,  PHASE , NDAT  A ,  LABEL .3.0,01 
550  J 1  ==  J  (  24  ) 

J2= J  <  25  > 

IF (MODL.NF  .  1  )  00  TO  565 
HO  563  I=J1,J2 
PHAZE=PHASL( I )*PI/180. 
sxi=sxi  +sn>KE:x(  i  > 

SX2-SX2+S I  OR LX ( I )**2 
St  1 =SY 1 T  PHAZE 
SXY  =SX Y  +  S I ORT  X ( I )*PHAZE 
563  NCUUNT  (3>=NCOUN7  (3>+l 

SLOPE® ABSf (SY1*SX1-SXY*N COUNT (3)  > /< SX1**2- SX2*NCUUNT  (3> >  • 

WRITE (8.910) SLOPE 

910  FORMAT'.  IX,*  SLOPE  =  *,  014.7,/) 

565  WRITE  <  4 ) NBA 1 A » 1 . /INCX 
DO  580  1  =  1  ,  NDAT A 

580  WRITE  (  4  ) F'1EL  IK  I  )  ,  STOREX <  I +  1 NCX  > 

WRITE (8, 920) 

920  FORMAT  <  IX. *CUMFT.ETED  CALCULATION  ANH  PLOT,  L  XTENUEB .  * ,  /  ) 

RETURN 

ENH 

C********* ***************************************** *********** 

SUHROUI  INE  HGRAPH(X«  Y,N,  ILi , NO ,NP , NS > 

DIMENSION  X(  1  >  ,  Y(  1  >  ,  IIi<25>  *  IF  (N0.L0.2)  CALL  PLOT  (  - 1 . 05  •  2 . 1  0  • i 
IF  (N(J.t.0.2)  GO  TO  30  S  IF  (NO.LI.O)  GO  TO  10 

CALL  SCALE (X, 7. ,N,1 )  $  CALI  SCALE ( Y , 5 . • N , 1 > 

10  CALI  PLOT  (0.  .  11  .  .2)  *  CALL  PI  01  ( 8 . 5 , 1  1 .  , 2. > 

CALL  PLOT  (8.5.0.  ,2)  *  CAl.l.  PLOT  (0.  ,0.  «2> 

CALI  PLOT  (1.35,1  .35,-3)  *  CALL  Pl. OT ( 0 .  ,11 . 30 ,  -2 ) 

IF  <  IIP  1>  .L0.9V9)  GO  TO  25 

CALL  PLOT < . 1 ,-. 1 ,  3)  t  CALL  PLOT ( 0 . • -2 . , -2 > 

CALI.  SYMBOL  (  .25,  .3,  .07,  UK  1  )  ,90.  ,20) 

CALL  SYMBOL ( .45, .3, .07,10(3) ,90. ,20) 

CALL  SYMBOL (  .6,5,  .3,  .07,  IIK  13)  ,90.  ,20) 

CALL  SYMBOL ( .85, .3, .07, IIK  15) ,90. ,20) 

CALL  SYMBOL.  ( 1 . 05 ,  .  3 ,  .  07  ,  l  IK  5 )  ,  90 .  •  20) 

CALL  SYMH(.IL<  1  .  15,  .3,  .07,  UK  7)  ,90.  ,20) 

CALI.  PLOT  (0.  ,0.  ,3)  $  CALL  PLOT  <  1 .25,0.  ,2) 

CALL  PLOT  <1  .25,2.  ,  2. )  *  CALL  PLOT  (  0  .  ,  2  .  •  -  2 ) 

CALL  PLOT (- . 1 , . t ,-3) 

25  CALL  PLOT . . . 

CALL  PLOT (0. ,-H. 30,-2)  *  CALL  PLOT ( -5 . 8 « 0 . • -  2 > 

CALL  Pl.01  (5.3,  .75,-3) 

CALL  AX  I S(0.  ,0. , I  IK  9 ) ,-20,7. , 90 . , X < NT  1 ) , X < N+7 ) ) 

IF  (NO. EC).  3)  GO  TO  27 

CAl.l  AXJS(0.  ,0.  »ILK  1 1 )  ,20,5.  ,100.  ,  Y(HI  1  )  ,  Y(Nf  2)  ) 

GCI  TO  30 

27  CALL  AXIS(0. ,0. , IIK  17) ,20,5. ,180. , Y (Nt 1 ) , Y(NI 2)  ) 

30  Y  <  N  I  2  )  =— Y  (  N  I  2  )  *  CALL  l.  1NE  (  Y  ,  X  ,  N,  1  •  NP  »N!»  > 

Y  ( N  I  2  >  --Y  (  N+2  >  $  CALL  Pl  OT  (  1  .85 . -2 .  .10  ,  -3  ) 

RETURN  *  I  NO 

SUBROLIT  INF  AXIS(XO,  YO»L,NC,Rl  .  ANG,  RM  I N  ,  DR  ' 

HI  MENS  I  ON  1.(1)  t  A  =  ANG*3. 14159/180.  *  I.IX== .  1  *l'.US  (  A )  1  HY== .  1*81  N  l  A  ) 
JC=ISIGN<  1 ,  Nl. )  *  NNC=I  ABS  ( NC )  1  R=.l  1  N=  1  t  X-XO  1  Y-YOi 
10  CALL  FLO  I  ( X. ,  Y ,  3  )  *  X==X+HX  *  Y-YIHY  %  CALL  PI01(X,Y,2) 

CAL  L  F'l  07  (X -.21  *BY*IC»Y+  . 2 1. *HX* I C ,  2  ) 

IF(M.L(>.5>  CALL.  PIOT(X-  .  42*H  Y*  TC  ,  Y  (  .  42*UX*IC  •  2  ) 

IF  (N. HO.  10)  CALI  PLOT  (  X-  .  70*I'Y*  TC « Y  t  .  70*HX*IC  •  2  ' 
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SB  40  = 
5050= 
5060= 

50/0 

5PB0= 
Ski  90  = 
5900  = 
S9 1  0  = 
5920  = 

5930 = 

59  10 
5950 
5960  = 
59/0= 
5980- 
5990= 
6000  = 
6010  = 
6020= 
6030  = 
6040  = 
6050 
6060= 

60  70= 
6000  = 
6090= 
6100- 
6110  = 
6120= 
6130= 
6140  = 
6150  = 
6160 
6170 
6180 
6190  = 
6200= 
6210= 
6220  = 
6230= 
6240  = 
6250= 
6260  = 
6270= 
6280= 
6290  = 
6300= 
6310  = 
6320= 
6330  = 
6340= 
6350= 
6360= 
6370= 
6380 
6390  = 
6400  = 
6410  = 
6420= 
6430  = 
6440- 
6450= 
6460  = 
6470  = 


N-M01KN.  10)41  S  R=R4 .  I  *  IF  <  R  .  t.  T  .  hi.  >  till  Til  .10 
A=ANG- < IC4 1 >*45.  *  DX=10.*DX  t  DY=10.*DY 
C-  .1  754 . 1  25*18  *  I.i=  .194. 35*1  C 

x=xo4L*ny-n*riY  *  y=yo4c*hy  t  n*nx 

R=AMAX1  ( ADS(F:MIN>  .  APS  <  RM I N4 DR*RL  >  )  1  l==Al  tlli  1  01R ) 

IR=  INK  ABS(R)  )  t  IF1R.LT.0.)  IR=-(IR41>  4  TF<-=  IR-MUD  (  IFc »  3  ) 

F\  1  -FCMIN/l  0 .  **  IFD  *  IIR1=DR/10.**IR  t  k=0. 

20  ENC0DF<7f 101 fS>R1  *  CALL  SYMBOL  <X»Y.  .07fSfAf7>  *  R1=RJ<DR1 
X  Xt-DX  t  Y  =  Y4IiY  t  R=R4 1  .  %  IF  (K  .LF'.KL  >  lid  1(1  20 
R=(RL- .  l*NNC>/2.  *  C=  .14. 5#  1C 

X  =  X04R*ItX-C*DY  *  Y  =  YO  F-F.*DY 4C*tiX 

CALL  SYMBOL  (XiY.  .1  pLi  ANGf  NNC  )  4  IF1IK.EO.O)  RETURN 
ENCOIILIS- 102.S)  4  CALL  SYMBOL  (999.  .999.  r  .  10.S. ANli.5) 

CALL  WMF  RL  <  X  .  Y.  A) 

F.NEOdF  <  3 . 1  03  f  S  >  Ik  *  CALL  SYMBOL < X f Y r . 07 f Sf ANG f 3 > 

101  FORMAT <F 7.2) 

102  F  OK'MAT  (5H  *10) 

103  FORMAT (13) 

RETURN  t  END 

C ********************************************** ******************* 

c: 

C  SUBROUTINE  SCAL.F  (HAT  A  .  LENGTH  .  N  .  K  ) 

C 

C  IDEAL  BATA  =  NT  2  DIMENSIONED  Afv.  Y  OF  DA  I A  TO  HE  SCALED 

C  INTFLjFR  N  -  NUMBER  OF  DATA  F'OINIS 

C  REAL  LENGTH  =  LENGTH  OF  THE  F'l.OT  AXIS  ( E  .  G .  IN  INCHES) 

C  INTEGER  K  =  UNUSED  PARAMETER  INCLUDED  F  OR  COMf'AT  I D U.  I  (  Y 

C  WITH  THE  EQUIVALENT  CALCOMF’  SUBROUTINE 

C 

C.  TFIE  FOLLOWING  VALUES  ARE  RETURNED! 

t 

1.  DATA1N41  )  =  ADJUSTED  DATA  MINT  MUM 

L  DAT  A  <  N42  >  =  ‘NICE*  SCALE  F  ACT  OF?  IN  DATA  UNITS 

C  PER  LENGTH  UNIT  <E.G.  VOL  TS/INCH ) 

C 

c ************ ***************************************v*Y*******«** 
SI.IDROUT  INF  SCALE < DATA f L ENG THf N f K> 

RIAL  DATA(N)f  LENGTH f  SE ( S ) 

DATA  SE  / 1 , »  2..  2 . 5 f  5.f  10.  / 

C  COMPUTE  THE  RAW  SCALE  FACTOR 

DM  IN==IIMAX  =  DAT  A  <  1  ) 

DO  10  1=1 fN 

IF  (DATA! I)  .LT.  DMIN)  DMIN  =  DA  I A ( 1  ) 

IF  <  DA  I A ( I > •  GT .  DMAX )  DMAX  =  DAIA(l) 

10  CONTINUE 

C 

C  EXCLUDE  TRIVIAL  ERROR  CASES 

C 

DATA <  N4 1  )  =  DMIN 
DATA! N42)  =  1.0 

IF  (LENGTH  ,LE.  0.0  .OR.  DMAX  . EQ .  DMIN  )  RETURN 
RAWSF  =  (UMAX  -  DMIN)  /  LENGTH 

=C  RAWSF  -  SF  MAN  I  *  10.  **  SFEXF'f  WHERE  1  .LE.  Sf-MANT  .IT.  10 

=  SF'EXP  =  AIN  1  (  AL0G101  RAWSF  )  ) 

=  IF  (  RAWSF  .LI.  1.0  )  SFEXF'  =  SF  EXP  -  1.0 

=  SF  MAN  I  =  RAWSF  *  10.0  **  <  SF  EXF‘> 

=C  LOCATE  NEXT  LARGER  ’NICE*  SCAl E  FACTOR 

DU  20  1=1 f 5 
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6480=20  IF  <  SF  ( I  )  .GT.  SF MART  )  GCI  FO  30 

6490-  F'RIN  I  *  *  *  SCALE!  SCAI  F  FACTOR  EKR’CIK  ...  ■  »  RETURN 

6500-30  SFrN  1  OR  =  SF(T)  *  10.0  »*  SF  l  XF' 

6510= 

6520=C  COMFUFE  All  JUS  r  E  n  uata  minimum 
6530- 

6540=  AMJMIN  --  A 1 N 7  (  UMIN  /  SFNICF  )  *  SF NICE 

6550=  IF  <  AIiJHIN  .01.  FiMIN  >  AIUMIN  -  AFUMIN  -  SFNICF 

6560=  IF  <  (UMAX  -  AIUMIN)  /  SFNICF  .  l.tl .  LENGTH  )  HO  TO  40 

6570  = 

658o=C  NFFIIi  TO  USE  I  HE  NEXT  LARGER  SCALE  TAG  FUR 

6590= 

6600=  IF  <  I  .LT.  5  )  SEN  ICE  =  SF(I41)  *  10.0  **  SI  I  XF' 

6610=  IF  <  I  .  EG .  5  )  SF  NICE  =  20.0  *  10.0  **  GFF.XF' 

66:>0=  AU.IMIN  =  AIN  I  (  UMIN  /  SFNFCFI  >  »  SF  NICE 

6630=  IF  (  AFUMIN  .01.  DMIN)  AFUM1N  =  AFUMIN  -  SENILE 

6640=40  CtINFINUE 

6650=  riATA(Ntl)  =  AFUMIN 

6660=  FiA  TA  (  N4  2  )  =  SEN  ICE 

6670=  Iv-L  TURN 

661)0=  END 

6690--C#» *****************  ******** ************************ *********** 
6700=  COMF1. EX  FUNCTION  CERF(7Z) 

6710=  COMI  'L.F  X  22  .  2 »  A  ,  A 1  »  A2  .  U  .  HI  »  U2  .  T  .  F  1 

6730=  2  =  22. 

6730=  IF  <  CAPS <  2  )  .C'E  .3 . 1G0T030 

6740=  J-0 . 

6750=  A =2 

6760=  U=? 

6770=  10  J=  J4 1 

£,700=  U=-/.*Z*CMF'LX<FLOAJ  <2*J-1  )  .0.  >*U 

6790=  F<=U  /CMF'EX  <  F  I  (I  A  I  (.1  )  .0.  )/CMFLX<Fl  (JAT  <2*. 14  1  >  .0,  ) 

6800=  A=AtU 

6810=  T F  < ,1 . HE  .  1  000  )  GO!  CI50 

6820=  I F ( CAHS ( U/A )  . GE  .  <  1  . E - 10 )  )  GFJ  TCI  10 

6830=  CEF\T=  <1  . 1  28379167 * 0 .  )*A 

6840=  FcFTUKN 

6850=  30  I F'  <  FcEAl .  <  22  )  .1.1 . 0 .  )  2  =  -2Z 

6860=  A2=(l.»0.) 

"’70=  I<2=2 

0-  F2=A2/U2 

6FJ90=  A  1  =  7. 

6900=  Eil=Z*24<0.5.0.) 

6910=  F  1=A1/U1 

6920=  J-l 

6930=  40  J= JF 1 

6«40=  A  Z*  A 1  4  CMF'l  X  (  FLOAT  <  J )  72  .  .0.  )*A2 

6950=  B  =  Z*61 FCMPLX (F  LOAT <  J) 72. -O. )*U2 

6960=  F  =  A  /  F< 

6970=  IF  ( J.G7  .1000) GOTO  50 

6980=  IF  <CAF(S<  (F  -FI  )/F)  .LT  .  <  1  .E~10>  ) (.1(11060 

6990=  A2=A 1 

7000=  U2=F.<1 

7010=  A 1 = A 

7030=  !<1=U 

»  7030=  F 1 =F 

7040=  Gfir()40 

7050=  50  UR  I IE <R. 99) 

7060=  99  FORMAK  •  ERROR  FUNCTION  ROUTINE  FUD  NUT  CONUFRGE  *> 

7070=  IER=1 

7080  RETURN 

7090-  60  F 1 =(0.5.0. )*CFXF< -2*7) *F 

7100=  CF  Rl  =t . 1203791 67*F1 

7110=  CERF'=1  .-CFIFi’l 
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7  J  ?0 
7130= 
7140  = 
7150  = 
7160  = 
7170  = 
7180  = 
7190  = 
7200  = 
10  = 
?0= 


I F  < REAL <  2Z )  . L.T .0.  )  CERF  =  -CERF 
70  RETURN 
E  NO 

C#*#*##  ********************************  *****************•*'***** 
COMF'LE'X  FUNCTION  f  RESL  <  X  ) 

COMPLEX  EYE, Z. CEKF 
E  YE=  <  0 . • 1 . )  *  F1=2.*ASIN< 1 . ) 

Z=SURT (F  T  >*X*< 1 .-EYE)/?. 

FRESt =<1 . +EYE ) /? . *CERF  < Z ) 

FRESL=CONJG( FREED 
RETURN  *  ENLi 
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